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v'was  often  found  to  remain  remarkably  stable  for  hours.  Significant  changes  in  the  size  distribution 
were  often  associated  with  changes  in  air  masses  and  meteorological  conditions.  Some  of  the 
changes  can  be  explained  qualitatively.  Others,  such  as  a  pronounced  double  peak  which  occurred 
occasionally  in  the  fine-particle  range,  are  difficult  to  understand. 

Measurements  of  the  activation  parameter  for  submicrometer  particles  were  also  made  and 
yield  values  in  the  range  of  0.1  to  O.SS. 

Using  the  measured  size  distribution,  electromagnetic  extinction  was  calculated  for  the 
wavelength  range  0.3  to  12  micrometers.  The  calculations  compare  favorably  to  the  measured 
scattering  coefficient  at  0.S5  micrometers  in  wavelength.  Particles  in  the  radius  range  0.1  to  0.3 
micrometer  made  the  largest  contribution  to  the  total  extinction  at  optical  wavelengths  on  nearly  all 
occasions. 

The  instrumentation  and  methods  which  were  developed  at  NRL  are  described  in  some  detail 
in  the  report. 
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PREFACE 

This  report  is  written  for  two  audiences,  with  the  first  being  our  sponsors  at  NRL,  ONR,  and 
NAVAIR,  whom  we  wish  to  provide  with  a  self-contained  document  of  our  accomplishments  in  the 
area  of  aerosol  measurements  over  the  last  several  years  without  their  need  to  consult  numerous  other 
publications  in  scientific  journals.  Sufficient  background  is  given  so  that  our  sponsors  can  place  our 
work  in  perspective  and  distinguish  NRL's  unique  contributions.  The  second  audience  consists  of  spe¬ 
cialists  in  atmospheric  aerosols.  The  third  section  of  the  report  is  a  review  of  the  instruments  and 
methods  developed  at  NRL.  Much  of  this  material  can  be  found  in  the  publications  referenced.  The 
latter  part  of  the  report  contains  the  results  obtained  on  three  field  deployments,  which  are  for  the  most 
part  not  yet  reported  in  the  open  literature.  A  reader  not  interested  in  the  measurement  techniques  or 
already  familiar  with  NRL’s  capabilities  may  want  to  skip  directly  to  the  results  and  the  discussion. 
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MEASUREMENTS  OF  ATMOSPHERIC  AEROSOLS: 
EXPERIMENTAL  METHODS  AND  RESULTS 
OF  MEASUREMENTS  OFF  THE  EAST  COAST 
OF  THE  UNITED  STATES 


INTRODUCTION 

Atmospheric  aerosols  are  the  population  of  particles  suspended  in  the  atmosphere  and  range  in 
radius  between  atmospheric  ions  (about  1  nm)  to  precipitation-size  particles  (about  1  mm).  Few  stu¬ 
dies  deal  with  this  entire  size  range,  and  the  use  of  the  term  "aerosol"  can  often  be  confusing  until  the 
subset  under  discussion  is  further  specified.  Furthermore  size  is  only  one  way  to  classify  aerosols.  For 
example,  a  meteorologist  may  refer  to  the  total  number  of  particles,  regardless  of  size,  as  condensation 
nuclei  (or  Aitken  nuclei),  or  the  meteorologist  may  refer  to  the  smaller  subset  which  is  active  in 
natural  clouds  and  fogs  at  low  supersaturations  (<1%).  The  particles  in  this  cloud-active  subset  are 
called  cloud  condensation  nuclei  (CCN).  Aerosols  can  also  be  classified  by  their  mode  of  formation, 
chemical  composition,  or  method  of  measurement.  All  of  the  properties  of  aerosols  are  important  in 
the  total  evaluation  of  aerosol  effects,  and  many  of  the  properties  are  interrelated.  Thus,  the  scattering 
and  absorption  of  electromagnetic  (EM)  radiation  by  aerosol  particles  depends  on  both  their  size  and 
composition,  as  does  their  nucleating  effectiveness  in  clouds  and  fogs. 

Some  reasons  aerosols  are  important  to  the  Navy  are  the  following: 

•  Extinction  of  EM  radiation  used  in  surveillance,  guidance,  and  communication  systems  is 
due  both  to  the  absorption  bands  of  gaseous  components  of  the  atmosphere  and  to 
scattering  and  absorption  by  aerosols.  Gaseous  absorption  is  more  predictable  and  better 
understood  than  scattering  and  absorption  by  aerosols  in  the  atmosphere.  Further 
research  on  the  geographical  distribution  of  aerosol  properties  is  necessary  in  order  to 
better  design  and  evaluate  systems.  A  capability  to  predict  aerosol  effects  during  military 
operations  is  also  important. 

•  Weather  modification  and  fog-visibility  prediction  depend  on  aerosols.  All  cloud  and  fog 
droplets  form  on  aerosol  particles;  thus  cloud  microstructure,  colloidal  stability,  and  visi¬ 
bility  are  largely  determined  by  the  aerosol  load  prior  to  cloud  formation.  Cloud  seeding 
is  an  effort  to  change  the  cloud  properties  by  changing  the  aerosol  on  which  it  forms. 

•  Artificial  aerosols  generated  by  pyrotechnic  devices  to  obscure  naval  operations  are  con¬ 
stantly  being  reevaluated  and  updated. 

•  Aerosols  affect  the  solar  radiation  balance  directly  by  scattering  and  absorption  and 
indirectly  by  changing  the  cloud  cover  and  the  reflectivity  (albedo).  Although  a  change  in 
the  radiation  balance  may  have  no  immediate  effects  on  military  operations,  it  is 
significant  in  projecting  global  climate  changes  and  could  affect  military  planning. 

BACKGROUND 

In  1977  the  field  of  aerosol  research  was  surveyed  to  determine  which  area  of  atmospheric  aero¬ 
sols  was  least  understood,  which  area  is  the  most  important  to  the  Navy,  and  where  the  NRL  research 
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effort  could  make  the  most  significant  contribution.  It  was  readily  concluded  that  little  was  known 
about  particles  smaller  than  0.2  um  in  radius.  Larger  particles  could  be  sized  with  optical  particle 
counters  and  impactors,  and  single  particles  bigger  than  I  /um  m  radius  could  be  chemically  analyzed 
using  special  techniques.  The  size  distribution,  chemical  nature,  geographical  distribution,  and 
mechanism  of  formation  were  all  poorly  understood  for  particles  smaller  than  about  0.2  um.  This 
understanding  was  particularly  poor  for  the  marine  environment  but  less  poor  for  the  urban  environ¬ 
ment,  where  small  particles  exist  in  much  larger  concentrations  and  where  a  larger  effort  in  pollution 
studies  has  been  supported  (Hidy,  1972,  and  others).  After  reviewing  the  few  submicromeier-size  dis¬ 
tribution  measurements  made  in  the  marine  environment,  Jaenicke  (1978)  stated,  "It  is  difficuli  to 
decide  whether  observed  differences  are  due  to  differences  in  time  and  location  of  the  measurements  or 
whether  the  distributions  are  biased  by  the  method  of  measurement." 

Prior  to  1976  NRL  had  extensively  studied  those  atmospheric  particles  (CCN)  which  are  active  in 
cloud  and  fog  formation  where  the  supersaturations  is  rarely  greater  than  1%.  During  the  early  1970s 
NRL  made  significant  advances  in  instrumentation  to  measure  CCN  and  made  measurements  of  CCN 
concentrations  and  supersaturation  spectra  at  many  geographical  locations  from  NRL  aircraft  and  ships 
From  these  data  much  information  about  the  geographical  and  altitude  distribution  of  CCN  was 
obtained.  Although  the  sizes  of  CCN  were  not  determined  precisely,  it  was  shown  that  CCN  were  in 
the  radius  range  of  about  0.02  to  0.09  /xm  and  were  composed  of  a  large  fraction  of  soluble  material 
(Twomey,  1965,  and  Hoppel,  1979).  The  composition  was  originally  thought  to  be  sea  salt  However, 
volatility  measurements  made  by  NRL  indicated  the  material  was  much  too  volatile  to  be  sea  salt 
(Twomey,  1968,  1971,  Dinger,  et  al,,  1970,  and  Hoppel  et  al.,  1973)  and  was  most  likely  to  be  sub¬ 
stances  such  as  ammonium  sulfate  formed  by  gas-to-particle  conversion  within  the  atmosphere  itself 
(Hoppel  et  al.,  1973,  and  Hoppel,  1975). 

The  dependence  of  visibility  on  relative  humidity  has  been  of  considerable  importance  within  the 
Navy.  This  humidity  dependence  is  caused  by  the  uptake  of  water  by  soluble  aerosols  at  relative  humi¬ 
dities  above  50%.  NRL  developed  a  theoretical  model  of  the  growth  of  particles  with  increasing  relative 
humidity  using  measured  data  on  the  chemical  composition  of  larger  particles  (Fitzgerald,  1975,  1978). 
This  model  has  been  used  in  a  number  of  Navy  programs.  More  recently  NRL  has  developed  a  unique 
method  to  measure  the  solubility  (hygroscopicity)  of  submicrometer  aerosols  (Hoppel,  1979,  and 
Fitzgerald  and  Hoppel,  1982).  This  method  will  be  discussed  later. 

With  NRL’s  historical  interest  in  CCN,  which  are  a  subset  of  the  submicrometer  particles,  it  was 
quite  natural  to  turn  our  attention  to  the  entire  population  in  that  difficult  and  little  explored  region 
below  0.3  /xm.  It  is  a  difficult  region  experimentally  because  the  particles  are  below  the  optically  detect¬ 
able  range  and  are  too  volatile  to  be  probed  by  conventional  methods  such  as  electron  microscopy. 
Since  the  particle  sizes  are  of  the  order  of  the  molecular  mean  free  path  (Fig.  1),  neither  continuum 
mechanics  nor  kinetic  theory  is  applicable,  and  special  theoretical  methods  are  required. 

In  1977  the  NRL  mobility  analyzer  was  developed  (Hoppel,  1978).  This  instrument  (to  be 
described  later)  has  given  us  a  unique  capability  to  size  particles  in  the  range  0.006  to  0.5  /ixm.  It  has 
been  used  in  two  ways:  primarily  as  a  size  spectrometer  and  secondarily  as  a  means  to  segregate  parti¬ 
cles  of  nearly  monodisperse  size  for  further  study.  It  was  first  used  on  a  USNS  Haye.7  erwse  in  1977 
but  with  limited  sensitivity;  even  so,  the  results  were  unique  (Hoppel,  1979).  In  1978  a  new  instru¬ 
ment  was  developed  (Hoppel  et  al.,  1980)  to  detect  and  count  single  particles  transmitted  by  the 
analyzer.  When  used  with  the  mobility  analyzer,  this  instrument  gave  us  an  extremely  sensitive  system 
which  could  be  automated.  The  results  of  size  distribution  measurements  to  be  given  in  this  report  are 
those  obtained  with  the  automated  system.  For  particles  larger  than  about  0.4  /xm  in  radius,  commer¬ 
cially  available  optical  particle  counters  were  used. 
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Fig.  1  —  Typical  si/e  distributions  of  condensation  nuclei 
(CN)  for  urban,  continental,  and  coastal  marine  environ¬ 
ments.  Also  shown  are  some  benchmark  dimensions 


Figure  1  illustrates  the  large  concentration  of  particles  in  the  size  range  below  0.2  as  opposed 
to  the  larger  size  range  accessible  to  optical  probes.  These  size  distributions  were  taken  with  the  NRL 
mobility  analyzer  and  are  typical  of  the  size  distributions  to  be  discussed  in  this  report. 

The  ability  to  segregate  small  particles  of  nearly  monodisperse  size  from  a  polydisperse  sample  is 
of  prime  importance  in  determining  particle  composition,  if  the  chemical  composition  and  dry  size  of  a 
particle  are  known,  then  according  to  Kohler’s  theory  it  is  possible  to  predict  its  size  at  higher  relative 
humidities  as  well  as  the  critical  supersaturation  at  which  the  particle  will  be  nucleated.  Simultaneous 
data  on  the  size  distribution  and  supersaturation  spectrum  of  known  pure-salt  aerosols  taken  at  the 
Second  International  Workshop  on  Cloud  and  Ice  Nuclei  in  1970  (Grant,  1971)  implied  a  relationship 
between  dry  size  and  critical  supersaturation  which  deviated  significantly  from  the  relationship  predicted 
by  Kohler’s  theory.  The  acceptance  of  Kohler’s  theory  and  its  extensive  use  in  cloud  physics  was  sub¬ 
sequently  questioned  in  the  literature  (Katz  and  Kocmond,  1973).  One  of  the  first  uses  of  the  mobility 
analyzer  was  to  segregate  a  narrow  size  range  of  particles  and  transmit  it  to  the  thermal-gradient 
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diffusion  cloud  chamber,  where  the  critical  supersaturation  for  nucleation  could  be  measured.  The 
results  agreed  with  those  predicted  by  Kohler's  theory  for  known  salts  (Gerber  et  al.,  1977)  and  were 
readily  accepted  as  validating  that  theory.  The  difficulty  in  the  earlier  experiments  was  that  the  meas¬ 
urements  were  much  less  accurate  than  was  claimed.  Extension  of  these  measurements  to  natural  aero¬ 
sols  has  led  to  the  first  measurements  of  the  solubility  of  natural  particles  smaller  than  0.1  pcm.  Some 
of  these  data  have  recently  been  published  (Fitzgerald  et  al.,  1982)  and  will  also  appear  in  this  report. 

In  1979  and  1980  a  number  of  field  measurements  were  made  with  the  automated  system  together 
with  supporting  measurements.  During  June  1979  12  days  of  measurements  were  made  on  the  Atlantic 
coast  at  NASA’s  Wallops  Island  rocket  launching  site.  The  purpose  of  these  measurements  was  to 
differentiate  the  size  distribution  of  marine  and  continental  aerosols.  In  July  1979  several  weeks  were 
spent  making  measurements  at  NRL's  Potomac  River  waterfront  to  test  the  system  and  to  sample 
urban  air.  Then  in  April  1980  we  participated  in  a  2-week  cruise  off  the  coast  of  the  United  States  with 
scientists  from  the  NRL  Environmental  Biology  Branch.  The  object  of  our  participation  was  to  study 
the  change  in  the  aerosol  size  distribution  in  air  masses  advecting  off  the  east  coast.  Then  in  the  fall  of 
1980  we  participated  in  the  International  Cloud  Condensation  Nuclei  Workshop.  In  addition  to  taking 
the  thermal-gradient  diffusion  cloud  chamber  (TGDCC)  and  the  newly  developed  isothermal  haze 
chamber  (IHC)  for  comparison  with  other  such  instruments,  NRL  was  invited  to  bring  the  NRL  mobil¬ 
ity  analyzer  to  characterize  the  size  distribution  of  aerosols  used  at  the  workshop. 

In  this  report  the  results  of  the  first  three  field  trips  mentioned  will  be  presented  with  emphasis  on 
the  data  obtained  on  the  1980  Hayes  cruise.  The  results  of  measurements  at  the  International 
Workshop  have  been  adequately  reported  in  the  workshop  proceedings  and  in  four  articles  by  NRL  per¬ 
sonnel  (Fitzgerald  et  al.,  1981,  Hoppel,  1981a,  1981b,  and  HoppcI  and  Wojchiechowski,  1981)  in  the 
special  issue  of  the  Journal  de  Recherches  Atmospheriques  devoted  to  the  workshop. 

The  primary  purpose  of  this  field  work  has  been  to  obtain  a  data  base  which  will  lead  to  an  under¬ 
standing  of  submicrometer  aerosols,  their  mode  of  formation,  and  the  mechanisms  which  shape  the  size 
distribution  and  are  responsible  for  the  dynamic  evolution  of  the  size  distribution  in  air  masses  advect¬ 
ing  over  the  ocean.  Emphasis  will  be  on  the  submicrometer  particles  as  opposed  to  the  la'-ger  sea-salt 
particles  which  have  been  the  subject  of  many  other  investigations.  New  data  on  the  solubility  of  sub¬ 
micrometer  particles  will  be  presented  as  will  supporting  data  such  as  visibility  and  radon  concentra¬ 
tions.  Extinction  of  EM  radiation  in  the  wavelength  range  0.3  to  12  /am  due  to  aerosols  is  calculated 
from  the  size  distributions  and  compared  to  the  scattering  coefficient  at  0.55  /am  measured  with  a 
nephelometer.  Instrumentation  which  has  been  developed  at  NRL  will  be  described  in  some  detail, 
since  it  represents  relatively  new  technology. 

INSTRUMENTATION  AND  METHODS 

The  goal  of  the  work  reported  herein  has  been  to  increase  our  understanding  of  atmospheric  aero¬ 
sols.  We  would  have  preferred  to  accomplish  this  goal  by  using  commercially  available  instruments. 
Unfortunately,  however,  one  reason  so  little  is  known  about  submicrometer  particles  has  been  the 
absence  of  a  measurement  technology.  During  the  past  few  years  much  effort  has  been  expended  in 
developing  a  capability  for  submicrometer  aerosol  measurement.  This  effort  has  been  profitable  and 
has  resulted  in  the  NRL  mobility  analyzer,  a  variable-humidity  nephelometer,  and  an  isothermal  haze 
chamber.  It  has  also  made  possible  the  unique  results  given  in  the  next  main  section  of  this  report  and 
in  recently  published  papers. 

The  Mobility  Analyzer 

The  size  distribution  of  small  particles  can  be  measured  with  an  electron  microscope  by  first  elec¬ 
trostatically  precipitating  the  particles  onto  an  clcctron-microscope  grid.  This  method  is  not  only  tedi¬ 
ous  but  is  of  doubtful  validity  wb,.  "t  to  atmospheric  aerosols  which  contain  a  large  fraction  of 
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volatile  material  that  will  not  survive  the  vacuum  environment  of  the  electron  microscope.  In  the 
atmosphere  most  size  distributions  of  particles  below  0.2  /am  have  been  obtained  by  measuring  the 
transmission  of  particles  through  a  diffusion  bank  as  a  function  of  airflow.  The  number  of  particles  lost 
to  the  walls  depends  on  the  diffusion  coefficient  of  the  particles.  The  inversion  of  the  transmission 
characteristic  gives  the  distribution  of  diffusion  coefficients,  from  which  the  size  distribution  can  be  cal¬ 
culated.  The  inversion  technique  essentially  involves  taking  two  derivatives  of  experimental  data,  a 
process  which  is  inherently  inaccurate.  Still  other  measurements  have  been  made  with  a  commercial 
instrument  called  the  electric  aerosol  analyzer,  originally  due  to  Whitby  (and  described  by,  for  example, 
Lundgen  et  al.,  1979).  This  instrument  sizes  particles  down  to  radii  below  0.01  um  and  has  been  used 
in  a  number  of  air-pollution  and  laboratory  studies  but  lacks  the  sensitivity  and  accuracy  required  in  our 
studies. 

The  NRL  mobility  analyzer  (Fig.  2)  is  a  system  for  measuring  the  size  distribution  of  particles 
between  0.006  and  0.5  ixm  in  radius.  The  aerosol  sample  is  first  brought  to  charge  equilibrium  by  pass¬ 
ing  the  sample  through  a  region  of  bipolar  ionization.  The  sample  air  enters  the  mobility  analyzer 
through  a  slit  in  the  outer  cylindrical  wall  and  is  confined  to  a  thin  laminar  layer  along  the  wall  by 
filtered  sheath  air.  A  small  amount  of  filtered  air  is  extracted  through  a  slit  in  the  inner  electrode. 
When  a  voltage  is  applied  to  the  inner  electrode,  those  particles  which  are  charged  and  which  lie  in  a 
narrow  mobility  range  will  be  withdrawn  with  the  air  extracted  through  the  inner  electrode.  The  parti¬ 
cles  in  the  extracted  air  are  counted  in  a  cloud  chamber  after  they  are  nucleated  and  have  grown  to 
optically  detectable  sizes.  From  these  data  on  the  particle  concentration  in  the  extracted  air  as  a  func¬ 
tion  of  voltage,  the  mobility  distribution  of  the  charged  fraction  is  obtained.  Knowing  the  mobility  dis¬ 
tribution  of  the  charged  fraction  (of  one  polarity)  and  knowing  that  the  particles  are  in  charge  equili¬ 
brium  and  have  a  Boltzmann  charge  distribution,  one  can  obtain  the  size  distribution  by  an  iterative 
procedure.  The  iterative  procedure  corrects  for  multiply  charged  particles  which  have  the  same  mobility 
as  smaller  singly  charged  particles.  The  theory  of  the  mobility  analyzer  and  the  analytical  method  for 
recovering  the  size  distribution  is  given  in  Hoppel  (1978). 


This  instrument  can  be  referred  to  as  a  differential  mobility  analyzer  of  second  order.  It  differs 
from  a  simple  aspirated  condenser  (for  example  a  Gerdien  chamber)  which  measures  integral  properties 
of  the  mobility  distribution  such  as  conductivity  and  total  ion  count.  In  the  case  of  an  integral  counter 
two  derivatives  of  the  experimental  cu'rent-voltage  curve  must  be  taken  to  obtain  the  mobility  distribu¬ 
tion  function.  Second-c  ;er  chamber.'  (such  as  the  NRL  mobility  analyzer)  are  chambers  where  both 
derivatives  are  taken  i:  ••’•■nentp’  >y  dividing  the  incoming  air  into  sheath  and  sample  air  and 
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measuring  the  ionic  current  to  only  a  small  segment  of  the  electrode.  In  a  second-order  chamber  the 
data  are  directly  related  to  the  mobility  distribution  itself. 

In  mobility  chambers  the  ions  are  usually  detected  by  the  electrical  current  generated  when  their 
charge  is  transferred  to  an  electrode  upon  contact.  For  atmospheric  aerosol  measurements  this  current 
is  too  small  to  be  measured;  therefore  second-order  chambers  have  not  been  used.  An  important 
feature  of  the  NRL  mobility  analyzer  is  that  the  electrical  current  is  not  measured  but  rather  the  parti¬ 
cles  are  removed,  nucleated  in  a  cloud  chamber,  and  counted  individually.  This  increases  the  sensi¬ 
tivity  vastly  over  that  obtainable  from  an  ionic  current  measurement.  (Some  history  on  this  type  of 
instrument  is  given  by  Knutson  and  Whitby,  1975,  and  Liu  and  Pui,  1974.) 

The  concentration  of  particles  in  the  extracted  air  is  much  lower  than  in  the  sample  air.  The  first 
data  from  the  mobility  analyzer  were  obtained  by  measuring  the  concentration  in  the  extracted  air  with 
a  Poliak  counter.  Although  the  long-tube  Poliak  counter  is  probably  the  most  sensitive  of  the  light¬ 
scattering  condensation-nucleus  (CN)  counters,  its  use  imposed  rather  severe  limitations  on  the  sensi¬ 
tivity  of  the  system.  Useful  size  distributions  were  obtained  only  when  particle  concentrations  were 
high.  When  the  total  count  dropped  below  about  1000  particles/cm\  the  system  could  not  be  used. 
Even  at  high  concentrations  the  sensitivity  was  not  adequate  to  detect  particles  in  all  1 1  of  the  size 
(mobility)  channels.  Despite  these  limitations,  this  early  version  demonstrated  the  usefulness  of  the 
mobility  analyzer  and  produced  useful  results  (Hoppel,  1979). 

In  an  attempt  to  increase  the  sensitivity  of  the  system,  a  single-particle  CN  counter  was  developed 
at  NRL.  This  CN  counter  now  allows  much  greater  sensitivity  and  has  made  it  possible  to  automate  the 
system.  Although  this  new  continuous  particle  counter  is  used  as  part  of  the  mobility  analyzer  system, 
it  is  novel  and  will  possibly  have  other  applications.  Therefore  the  new  continuous  CN  counter  will  be 
described  in  some  detail  here. 

Figure  3  shows  the  design  of  the  new  single-particle  CN  counter.  The  sample  air  enters  along  the 
axis  of  the  chamber  and  is  surrounded  by  filtered,  humidified  sheath  air.  The  walls  are  water  saturated; 
alternate  segments  are  maintained  at  alternating  temperatures.  The  basic  principle  is  similar  to  that 
employed  in  the  TGDCC,  where  a  supersaturation  exists  between  horizontally  oriented  hot  and  cold 
plates  with  saturated  walls.  The  TGDCC  is  limited  to  measuring  panicles  activated  in  a  narrow  range  of 
supersaturations  between  about  O.Z'Mi  and  1.0%  because  of  well-documented  reasons  related  to  fallout 
during  the  time  required  to  establish  the  equilibrium  supersaturation  in  the  TGDCC.  The  advantage  of 
the  segmented  chamber  is  that  particles  falling  vertically  are  not  lost  and  the  segmented  geometry 
suppresses  large-scale  convective  motions  which  would  occur  between  vertically  oriented  hot  and  cold 
walls.  The  criterion  in  the  segmented  chamber  is  that  the  airflow  rate  must  be  such  that  the  residence 
time  in  a  single  segment  is  short  compared  with  the  time  required  for  moisture  or  heat  to  diffuse  in 
from  the  walls  to  the  center,  whereas  the  residence  time  of  the  sample  in  the  entire  chamber  must  be 
longer  than  the  diffusion  time.  When  this  criterion  is  met.  the  air  along  the  axis  will  reach  a  tempera¬ 
ture  and  vapor  pressure  which  is  approximately  midway  between  the  temperature  and  vapor  pressure  of 
each  segment.  Along  the  axis  the  asymptotic  supersaturation  is  the  same  as  that  existing  in  the  center 
of  a  TGDCC  working  at  the  same  temperature  difference.  Figure  4  shows  the  development  of  the 
supersaturation  within  the  chamber.  Since  the  sheath  air  is  prehumidified,  so  that  the  supersaturation 
begins  at  0%,  the  upper  curves  are  applicable  to  our  mode  of  operation.  Details  of  the  design  and 
analysis  of  the  chamber  appear  in  Hoppel  et  al.  (1979). 

In  practice  the  temperature  difference  is  set  such  that  the  asymptotic  (maximum)  supersaturation 
is  about  20%,  so  that  all  atmospheric  particles  in  the  size  range  of  interest  are  nucleated  and  grow  larger 
than  1  |im  in  radius.  The  particles  are  then  individually  counted  by  a  commercially  available  optical 
particle  counter  (a  Royco  optical  particle  counter,  which  was  used  in  this  application  only  to  count  the 
activated  droplets  and  not  to  size  them). 
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DISTANCE  ALONG  THE  AXIS  Icml 

Fig  4  —  ncvcIopnK'nl  of  supersaturation  in  the  segmented  chamber  iFig  '3i  I  he  radial  distance  in  the 
chamber  is  indicated  as  the  normalized  distance  R  such  that  /?  =  n  is  along  the  avis  and  -  1  is  along  the 
wall. 


In  the  winter  of  1979  the  mobility  analyzer  was  automated  using  an  HP-9825  desktop  calculator 
and  appropriate  interfaces.  A  16-bit  binary  interface  was  used  to  control  a  high-voltage  power  supply  (0 
to  5000  V)  which  supplied  the  polarizing  voltage  to  the  mobility  analyzer.  Another  interface  controlled 
and  read  the  Royco  optical  particle  counter  at  the  output  of  the  segmented  chamber,  A  real-time  clock 
in  the  HP-9825  controlled  the  timing  such  that  the  voltage  cycled  through  12  voltage  steps  between  0 
and  5000  V.  After  each  voltage  change,  60  s  are  allowed  for  a  new  equilibrium  to  be  established,  and 
then  two  consecutive  12-s  counts  are  obtained  from  the  particle  counter  for  each  voltage.  The  average 
of  the  counts  is  stored,  and  at  the  end  of  the  run  the  size  distribution  is  calculated,  printed  out,  and 
stored  on  a  magnetic-tape  cartridge.  A  plotter  has  been  added  to  the  system  so  that  the  size  distribu¬ 
tion  can  be  plotted  after  each  run.  A  single  size  distribution  takes  about  15  minutes.  The  increased 
sensitivity  of  single-particle  counting  now  allows  particles  to  be  counted  in  all  1 1  channels  under  all 
conditions  encountered  to  date. 

In  field  work  we  have  consistently  used  only  II  channels  to  cover  the  stated  size  range.  More 
channels  could  be  used  to  improve  the  size  resolution,  but  this  would  increase  the  time  required  to 
obtain  the  size  distribution.  The  selection  of  1 1  channels  has  been  a  good  compromise  between  resolu¬ 
tion  of  sizes  and  resolution  of  temporal  changes. 

The  size  distribution  n(r)  is  defined  as  the  number  of  particles  dNin  a  radius  interval  dr. 


In  calculating  the  size  distribution  from  the  mobility  analyzer  data,  the  infinitesimals  of  difTcrential  cal¬ 
culus  must  be  replaced  by  finite  differences.  The  procedure  is  given  in  Hoppel  (1978).  The  size  distri¬ 
bution  is  usually  plotted  on  log-log  paper  and  can  be  displayed  in  a  number  of  alternate  ways  depending 
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on  the  application.  Since  nir)  decreases  so  rapidly  with  increasing  size  (Fig.  1),  it  has  become  quite 
common  to  plot  the  logarithmic  derivative  instead; 


dN 

rfdog  r) 


=  In^dO)  r  ~  =  2.303  r 

dr  dr 


(2) 


This  form  of  the  size  distribution  gives  more  weight  to  the  larger  particles  and  hence  makes  the  distri¬ 
bution  less  steep. 


The  area  distribution  is  intimately  related  to  optical  properties  of  aerosols  and  is  defined  as  the 
size  distribution  weighted  by  the  geometric  cross  section: 

=  IT  r^  nir).  (3) 

dr 

The  integral  of  the  area  distribution  function  over  all  radii  is  the  total  cross-sectional  area  presented  by 
the  entire  population  of  particles  in  1  cm^. 


The  manner  in  which  the  volume  (or  mass)  is  distributed  among  the  various  sizes  is  given  by  the 
volume  distribution 


dr 


Integration  of  the  volume  distribution  over  all  sizes  gives  the  total  volume  in  1  cm^. 


(4) 


In  the  area  and  volume  distributions  the  larger  particles  are  given  increasingly  more  weight  and 
hence  skew  the  distribution  toward  larger  sizes.  In  this  report  log-log  plots  of  size  distributions 
idNidr)  and  area  distributions  idA/dr)  will  be  the  display  forms  used. 

Thermal-Gradient  Diffusion  Cioud  Chamber  (TGDCC) 


The  TGDCC  is  an  instrument  which  records  the  concentrations  of  particles  nucleated  at  low 
supersaturations  (0.2  to  1%).  Since  supersaturations  in  natural  clouds  rarely  exceed  1%,  the  TGDCC 
measures  the  subset  of  particles  upon  which  natural  clouds  form.  This  subset  of  condensation  nuclei 
(CN)  is  referred  to  as  cloud  condensation  nuclei  (CCN).  The  TGDCC  was  developed  largely  at  NRL 
in  the  late  1960s  by  Twomey  (1967)  and  is  now  widely  used  to  measure  CCN.  In  addition  to  using  the 
instrument  to  obtain  the  supersaturation  spectrum  between  0.2'Ki  and  1.0''/(i  supersaturation,  the 
instrument  was  used  in  conjunction  with  the  mobility  analyzer  to  obtain  information  on  the  solubility  of 
small  particles  during  the  field  operations.  This  novel  use  of  the  instrument  will  be  described  in  the 
next  subsection.  Since  the  NRL  TGDCC  has  been  described  in  several  papers  (Hoppel  and 
Wojciechowski,  1981),  only  a  brief  description  is  given  here. 


The  TGDCC  chamber  consists  of  two  plates  7.5  cm  in  diameter  separated  by  1.25  cm  and  covered 
with  saturated  filter  paper.  A  cylindrical  glass  wall  at  the  circumference  of  the  plates  and  separating 
them  forms  the  chamber.  The  top  plate  is  at  room  temperature,  and  the  bottom  plate  is  cooled  with  a 
thermoelectric  cooler.  The  temperature  difference  is  measured  with  several  sets  of  thermocouples 
imbedded  in  the  plates.  A  high-intensity  mercury-arc  lamp  with  appropriate  lens  and  collimating  slits 
which  are  imaged  in  the  chamber  produce  an  illuminated  volume  which  is  viewed  at  a  90°  scattering 
angle  by  an  8-mm  camera  and  a  video  camera  system. 


The  video  system  permits  immediate  playback  and  stop-frame  counting  of  the  cloud  droplet  con¬ 
centration  in  the  chamber.  The  30-cm  (12-in.)  monitor  screen  is  marked  off  to  identify  the  viewing 
dimensions  for  ease  in  counting  the  droplet  images. 


Air  to  be  analyzed  is  drawn  into  a  10-liter  holding  chamber  by  drawing  down  a  metalized  Mylar 
diaphragm.  This  air  is  then  continuously  circulated  through  the  holding  chamber  until  a  measurement 
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is  made.  To  make  a  measurement,  the  eireulating  air  is  shut  off,  and  a  slight  pressure  is  applied  to  the 
diaphragm  to  foree  air  from  the  holding  chamber  to  the  TfiDCf  The  count  is  recorded  both  on  sideo- 
lape  and  a  number  of  frames  of  8-mm  mosie  film  I  he  8-mm  Bolex  camera  is  operated  m.inualU  while 
the  chamber  is  \iev\ed  through  the  camera  I  he  first  frame  is  taken  1  to  2  s  after  the  \al\e  is  closed, 
and  ten  to  15  frames  are  taken  manuall>  until  it  is  clear  that  the  maximum  particle  count  has  passed 
Two  or  three  successise  samples  are  recorded  at  each  of  four  supersaiurations  (as  established  b)  chang¬ 
ing  the  temperature  differences  between  the  plates) 

The  video-tape  count  is  quicker  and  easier  to  obtain  and  is  the  count  given  in  this  report  How¬ 
ever,  it  has  always  been  our  experience  with  atmospheric  aerosols  that  the  photographic  count  is  always 
about  10  to  15'yn  higher  We  attribute  this  to  the  smaller  detectable  si/e  obtained  with  the  photographic 
method.  The  minimum  detectable  si/e  for  our  photographic  system  is  estimated  to  be  about  0.5  /am, 
whereas  the  minimum  detectable  si/e  for  the  video  system  is  estimated  to  be  just  under  1  /am. 

There  are  definite  limits  on  the  supersaturation  range  for  which  the  results  of  the  TGDCC’  are 
valid.  The  range  of  validity  is  usually  given  as  0.2  to  l"/i.  However,  even  within  this  range  there  can 
be  significant  errors,  depending  on  the  nuclei  spectrum  (si/e  distribution)  being  measured.  The  accu¬ 
racy  limitations  of  TfiHC'Cs  m  general,  and  of  our  chamber  in  particular,  were  investigated  bv  Hoppel 
and  Wojciechowski  (1076) 

The  smaller  particles  which  have  critical  supersaturations  about  the  same  as  the  maximum  super¬ 
saturation  in  the  TCiOC'C  are  not  nucleated  until  the  chamber  has  reached  its  steady -state  value  and 
then  grow  more  slowly  than  the  larger  particles  which  are  nucleated  before  the  chamber  has  reached 
equilibrium.  The  object  is  to  find  an  interval  of  time  when  the  less  active  particles  have  grown  to 
minimum  detectable  si/e  and  the  larger  (more  active)  particles  have  not  yet  started  to  fall  out,  .X  time 
which  satisfies  both  criteria  may  not  exist  for  all  si/e  distributions.  If  such  an  interval  of  lime  does 
exist,  then  it  should  evidence  itself  by  a  plateau  in  the  curve  of  number  of  panicles  detected  versus 
time.  A  well-defined  plateau  was  not  found  to  exist  in  the  data  on  natural  continental  aerosols 
presen'ed  by  Hoppel  and  Wojehiechowski  (W76).  .Mofs  and  Carstens  (l'^76)  did  a  numerical  simula¬ 
tion  of  the  TCiOCC  which  predicted  large  errors,  depending  on  the  nuclei  distribution  and  minimum 
detectable  size. 

Our  method  for  obtaining  the  count  was  as  follows:  two  (or  three)  successive  samples  were 
recorded  at  each  supersaturation.  The  recordings  from  these  samples  were  played  back  and  visually 
examined  to  obtain  the  succession  of  frames  for  each  sample  where  the  maximum  count  occurred 
Several  of  these  frames  were  counted,  and  the  maximum  count  was  obtained.  The  maximum  counts 
from  two  (or  three)  samples  were  averaged  to  obtain  the  CCN  count  given  in  this  report. 

Measurements  of  the  Hygroscopicity  (Solubility)  of  Submicrometer  Particles 

Theory 

The  relationship  between  the  saturation  ratio  S  (relative  humidity/ 100)  of  the  air  and  the  equili¬ 
brium  radius  rof  an  aqueous  solution  droplet  may  be  expressed  as  (Pruppacher  and  Klett,  1978) 


where  is  the  water  activity  of  the  aqueous  solution,  rr  is  the  surface  tension  of  the  droplet.  /?,  is  the 
specific  gas  constant  of  water  vapor,  and  T  is  the  temperature.  The  water  activity  describes  the  reduc¬ 
tion  in  vapor  pressure  due  to  the  dissolved  material  and  is  equal  to  the  equilibrium  saturation  ratio  over 
a  plane  surface  of  the  solution.  The  exponential  factor  in  Kq.  (5)  describes  the  increase  in  vapor  pres¬ 
sure  above  that  of  a  plane  surface  due  to  the  curvature  of  the  droplet. 
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Atmospheric  aerosol  particles  consist  of  water-soluble  and  water-insoluble  substances,  the  soluble 
component  normally  being  a  complex  mixture  of  different  salts  (Winkler,  1973,  Junge  and  McLaren, 
1971,  and  Hanel,  1976).  For  droplets  formed  upon  such  particles  we  may  write  (Hanel,  1976) 

MJ m^)  /,x 

e  W'  s  m  ^ 

where  T-  is  the  mean  number  of  moles  of  ions  and  undissociated  molecules  per  mole  of  salt  mixture,  <t> 
is  the  observed  practical  osmotic  coefficient,  «  is  the  mass  fraction  of  soluble  material  in  the  aerosol 
particle,  is  the  mean  molecular  weight  of  the  salt  mixture,  is  the  molecular  weight  of  water,  wq 
is  the  dry  mass  of  the  particle,  is  the  mass  of  water  condensed  on  the  particle,  and  y  is  correction 
factor  taking  into  account  absorption  (or  desorption)  of  ions  by  the  insoluble  portion  of  the  particle, 
which  is  assumed  to  be  totally  submerged  in  the  aqueous  solution. 


We  denote  the  coefficient  of  mj in  Eq.  (6),  which  is  referred  to  as  the  exponential  mass- 
increase  coefficient  (Hanel,  1976),  by  tj: 


V  4)  f.  y 


(7) 


Thus  in  terms  of  17  the  expression  for  becomes 

-  n  (mJ ni^.) 


(8) 


The  coefficient  tj  is  a  function  of  solute  concentration  and  can  be  expressed  as  a  function  of  either 
or  mij m^.  If  we  make  the  simplifying  approximation  that  the  volume  of  the  solution  drop  is  equal  to 
the  sum  of  the  volumes  of  the  dry  particle  and  of  the  pure  water  condensed  on  it,  we  can  write 

=  1  +  (9) 


rl 


Ph-Wq 


where  tq  is  the  equivalent  volume  radius  of  the  dry  particle  and  po  and  Pn  are  the  densities  of  the  dry 
particle  and  water  respectively.  Using  Eq.  (9)  to  express  in  terms  of  r/ro  and  substituting  the 

result  into  Eq.  (8)  leads  to 


Ow 


(10) 


where  B  =  pov/Pw  =  Pof'<f>  *7  ^w/Pw^s-  (o  ^^s  the  aerosol  hygroscopicity  parameter.  The 

value  of  T),  and  hence  the  value  of  B,  is  not  constant,  owing  to  the  variation  of  4>  '''ith  the  concentra¬ 
tion  of  the  solution.  We  can  express  tj  and  5  as  functions  of  mJ or  a^. 


Since  tj  and  B  are  related  by  the  constant  factor  po/p».,  we  can  write 

(11) 

r? 

where  Tj°_and  are  the  values  of  tj  and  B  for  an  infinitely  dilute  solution,  that  is  for  «»,  =  1 .  For 

Oh,  -  1,  0=1. 

Hanel  and  Lehmann  (1981)  have  determined  (via  Eq.  (8))  tj  as  a  function  of  water  activity  for 
atmospheric  aerosol  samples  by  measuring  the  relative  mass  increase  (m^/nto)  as  a  function  of  o,  using 
a  microbalance  system.  Figure  5  is  based  on  their  measurements  and  shows  the  ratio  tj/tj"  as  a  function 
of  a*,  for  an  industrial  aerosol  sample  (curve  4)  and  several  samples  composed  of  background  aerosol 
combined  with  particles  from  industrial  and  local  sources.  The  key  to  the  curves  (the  list  in  the  figure 
title)  gives  the  values  of  tj"  and  in  addition  to  the  chemical  composition  of  the  soluble  portion  of 
each  of  the  samples.  To  obtain  the  value  of  tj",  Hanel  and  Lehmann  used  the  mixture  rule  for  semi¬ 
ideal  mixtures  of  electrolytes  to  extrapolate  the  value  of  tj  measured  at  =  0.96.  Several  features  of 
these  curves  are  noteworthy.  First,  the  curves  show  a  relative  maximum  in  tj/tj"  in  the  range 
0.75  ^  ^  0.8.  At  values  of  above  this  maximum  the  soluble  material  can  be  presumed  to  be 
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Fig  5  —  Normalized  exponential  mass-increase  coefficieni  rj/'l*’ 
versus  water  activity  a,  for  five  continental  aerosol  samples, 
based  on  measurements  by  Hilnel  and  Lehmann  (1981)  for 
increasing  relative  humidity.  The  chemical  compositions  of  the 
soluble  components  of  the  samples  along  with  the  t)°  ^nd  B® 
values,  keyed  to  the  curves,  are  as  follows: 

1.  0.57  (NH4)2  SO4  +  0,2J  NH4NOJ  -t  0.20  H2SO4, 

7)®- 0.12  and  fl®  -  0.23; 

2.  0.57  sea  salt  +  0.43  NH4  NO,; 

I,®  -  0.27  and  B®  -  0.90; 

3.  0  93  (NH4IJ  SO4  +  0.07  Hj  SO4; 

I)®  -  0.32  and  B®  -  0  59; 

4.  0.76  N^NO,,  +  0.24  H}S04; 

T,®  -  0.18  and  B®  -  0.34; 

5.  0.5  NH4NOJ  +  0  37  sea  salt  +  0  13  (NH4)j  SO4, 
n®  “  0  24  and  B®  -  0  61; 

6.  average  of  curves  t  through  5. 

totally  dissolved.  Second,  above  a»  =  0.78,  the  value  of  tj/t)*’  at  a  given  water  activity  shows  a  maxi¬ 
mum  variation  of  only  30®/o  and  a  deviation  of  ±15%  from  the  mean  curve.  The  value  of  on  the 
other  hand,  is  known  to  vary  by  as  much  as  about  a  factor  of  5,  due  to  variations  in  «,  po,  Wj,  and  f. 
Third,  for  a  given  sample  (a  given  mixture  of  electrolytes)  the  value  of  rj/T/"  is  fairly  constant  over  the 
range  0,78  <  <  0.99.  Thus  over  the  range  0.78  ^  a»  <  0.99  the  value  of  as  a  function  of  a» 

can  be  obtained  with  reasonable  accuracy  from  Eq.  (11)  by  measuring  the  value  of  and  using  the 
mean  curve  of  r)/r}°  versus  in  Fig.  5  or,  even  more  simply,  by  using  B  =  0.8  fl®. 

In  computing  particle  size  as  a  function  of  relative  humidity  from  Eqs.  (5)  and  (10).  the  depen¬ 
dence  of  surface  tension  on  the  concentration  of  the  solution  must  be  taken  into  account.  The  surface 
tension  of  an  aqueous  solution  can  be  described  by  (Hanel,  1976) 

nto 

u  =  <Th  +  - ,  (12) 

m.. 
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where  (r»  is  the  surface  tension  of  pure  water  and  b  is  constant  for  each  solute 
scopicity  parameter  B,  <r  may  be  expressed  as 

B  In 

<r  =  <r  «  ~  - r - , 

“  B 

where  B  =  bfpo/f)^.  The  surface  tension  of  water  can  be  approximated  by 

<rJT)  =  117.4  -  O  ISST, 
where  o-»  is  in  micronewtons  per  millimeter  (dynes  per  centimeter)  and  Tis  in  kelvins. 

An  error  in  a-  does  not  have  a  large  influence  on  the  equilibrium  size  of  a  particle  until  5  exceeds 
0.95  if  tq  >  0.1  or  until  S  exceeds  0.9  if  kq  is  as  small  as  0.05  /^m.  However,  at  these  high  values  of  S 
the  second  term  in  Eq.  (13)  is  sufficiently  small  compared  to  tr,  that  an  uncertainty  in  the  value  of  B 
of  as  much  as  50%  will  result  in  an  error  in  <t  of  less  than  5  to  10%.  Hanel  and  Lehmann  (1981)  have 
also  determined  the  values  of  i,  «,  and  po  for  their  atmospheric  aerosol  samples.  The  value  of  B  calcu¬ 
lated  from  these  data  ranges  from  9.4  to  47.7  pN/mm  and  has  the  average  value  of  20.2  pN/mm.  This 
value  of  B  >s  reasonable  one  to  use  in  computing  the  equilibrium  size  of  atmospheric  aerosol  particles  as 
a  function  of  relative  humidity. 


In  terms  of  the  hygro- 

(13) 


At  values  of  the  saturation  ratio  at  and  above  unity,  a  solution  droplet  may  be  assumed  to  be 
sufficiently  dilute  (such  that  a  ~  and  B  ~  fl®)  so  that  r^/rg  -  1  can  be  approximated  by  rVro  and 
the  exponential  terms  in  Eqs.  (5)  and  (10)  can  be  expanded  to  yield 


5  = 


r 


(14) 


where  A  =  2(t  J p„R^T  ~  4,33  x  10’’  trJT. 


The  relationship  between  the  dry  size  and  critical  supersaluration  (S,)  of  an  aerosol  particle  is 
obtained  by  differentiating  Eq.  (14)  and  setting  9S/dr  =  0,  yielding 


S,  = 


4/f^ 


21 B^ 


l'2 


'■o 


3  2  _ 


=  AVo 


3  2 


(15) 


Since  the  parameter  (=  potj°/Pw)  determines  the  supersaturation  at  which  particles  of  a  given  size 
are  activated,  it  will  be  referred  to  as  the  particle  activation  parameter.  In  this  study  Eq  (15)  was  used 
to  determine  the  values  of  S®  and  K  of  natural  aerosol  particles  from  measurements  of  the  critical 
supersaturation  and  dry  size  of  the  particles.  Knowing  the  value  of  fl",  we  computed  the  equilibrium 
size  of  the  aerosol  particles  as  a  function  of  relative  humidity  (over  the  range  0.8  ^  S  0.99)  from  Eqs. 
(5),  (10),  and  (ID,  using  the  approximate  values  B  =  0.8fi®  and  B  —  20.2  pN/mm. 


Experimental  Method 

The  critical  supersaturation  of  particles  of  known  size  was  measured  by  using  the  mobility 
analyzer  in  tandem  with  the  thermal-gradient  diffusion  cloud  chamber  (TGDCC).  The  mobility 
analyzer  is  designed  to  measure  the  size  distribution  of  aerosol  particles  and  is  described  elsewhere  in 
this  report.  However,  its  use  in  these  experiments  is  in  principle  the  same  as  the  electrical  aerosol 
classifier  described  by  Liu  and  Pui  (1974). 

The  first  step  in  the  experimental  procedure  is  to  bring  the  aerosol  sample  to  charge  equilibrium 
by  means  of  a  radioactive  source.  The  sample  is  then  passed  through  the  mobility  analyzer.  Particles 
in  a  preseDcted  narrow  size  range  are  segregated  from  the  total  aerosol  population  by  the  mobility 
analyzer  and  are  transmitted  to  the  TGDCC,  where  the  critical  supersaluration  of  the  particles  is  meas¬ 
ured.  This  method  was  employed  by  Gerber  et  al.  (1977)  to  verify  Kohler's  theory  for  pure  salt  parti¬ 
cles  of  known  composition. 
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To  determine  the  critical  supersaturation  of  particles  transmitted  by  the  mobility  analyzer,  the 
supersaturation  in  the  TGDCC  is  slowly  increased  by  gradually  raising  the  temperature  difference  AT 
between  the  plates  of  the  chamber.  Typically  five  to  ten  determinations  of  the  number  of  activated 
nuclei  are  made  in  the  time  it  takes  A  7"  to  increase  by  0.25°C.  The  average  number  of  particles 
activated  in  each  0.2S°C  interval  of  AT  is  then  plotted  as  a  function  of  the  mean  supersaturation  for 
that  temperature  interval.  Averaging  is  necessary  to  obtain  a  statistically  significant  count,  because  of 
the  small  number  of  particles  in  the  narrow  size  range  transmitted  to  the  TGDCC  and  because  of  the 
srhall  sensitive  volume  (0.067  cm^)  used  in  the  TGDCC.  Next,  the  number  of  multiply  charged  larger 
particles,  which  are  transmitted  along  with  the  more  numerous  singly  charged  particles  of  the  desired 
size,  is  estimated  from  the  particle  size  distribution  (which  is  measured  immediately  before  or  after 
each  measurement)  and  is  subtracted  from  the  number  activated  at  each  supersaturation.  The 
corrected  number  of  particles  at  each  supersaturation  is  then  normalized  to  the  maximum  number 
activated. 

As  an  illustration  of  the  method,  Fig.  6  shows  the  normalized  number  of  activated  nuclei  as  a 
function  of  supersaturation  for  a  run  made  on  30  July  1979.  The  transmitted  size  range  in  this  case  was 
from  0.0427  to  0.052  /^m  in  radius.  Particle  size  is  computed  from  particle  mobility  on  the  assumption 
that  the  particles  are  singly  charged.  The  range  of  critical  supersaturation  of  the  particles  in  the 
transmitted  size  interval  is  taken  to  be  between  the  points  where  10%  and  90%  of  the  particles  are 
activated.  These  points  are  identified  by  the  arrows  on  the  curves.  (Ninety  percent  rather  than  100% 
was  chosen  because  statistical  fluctuations  in  the  number  of  droplets  in  the  small  viewing  volume  used 
in  the  TGDCC  make  it  impossible  to  determine  the  exact  supersaturation  at  which  all  particles  are 
activated.)  Critical  supersaturation  is  plotted  against  particle  size  as  shown  by  the  blocks  in  Fig.  7  for 
the  three  transmitted  size  ranges  used  in  the  run  of  30  July.  The  widths  of  the  blocks  are  the  transmit¬ 
ted  size  intervals,  and  the  heights  of  the  blocks  are  measured  ranges  of  critical  supersaluration.  The 
values  of  critical  supersaturation  and  particle  radius  at  the  center  of  a  block  are  used  to  compute  the 
values  of  ff”  and  /f  from  Eq.  (15).  For  reference  Fig.  7  also  shows  the  theoretical  relationship  between 
critical  supersaturation  and  dry  radius  for  pure  ammonium  sulfate  particles  and  for  particles  character¬ 
ized  by  =  0.10. 
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Fig  6  —  Normalized  number  of  activated 
nuclei  as  a  function  of  supersaturation  for 
the  experimental  run  made  on  30  July  1979 
(7  -  31“C  1 
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continuous  lines  are  theoretical  curves  lor  ammonium 
sulfate  particles  and  for  particles  characterised  by  3 
fl“-0  10.  “ 


The  results  of  the  solubility  measurements  (B®  measurements)  will  be  discussed  later  in  this 
report. 


Calculation  of  Aerosol  Extinction  in  the  Wavelength  Range  0.3  to  12  /urn 

Electromagnetic  extinction  at  any  wavelength  depends  on  the  aerosol  size  distribution,  the  nature 
of  the  aerosol,  and  the  relative  humidity.  The  relative  extinction  at  different  wavelengths  largely 
depends  on  the  shape  of  the  aerosol  size  distribution  and  the  variation  of  the  complex  refractive  index 
with  wavelength. 

The  general  expre.ssion  for  the  aerosol  extinction  coefficient  is 

kf  =  irr^  Qeim,  a)  n(r)  dr,  (16) 

where  r  is  the  particle  radius,  n(r)  =  dN/dr  is  the  number  of  particles  per  cubic  centimeter  per 
micrometer  radius  interval,  and  is  the  Mie-theory  extinction  efficiency  factor.  Q,  is  a  function  of 
the  particle  size  parameter  a  =  lirrlx  and  the  complex  index  of  refraction  of  the  particle  m  —  n  —  ik. 
Exact  values  of  0,  are  computed  from  Mie-theory  codes.  However,  to  facilitate  the  calculation  of  kf, 
Fitzgerald  (1979)  developed  approximation  formulas  for  as  a  function  of  normalized  size  parameter 
p  =  2a(n  -  1).  For  aerosols  obeying  a  Junge  size  distribution  of  the  form  n(r)  =  Ar  ' ,  the  values  of 
k,  calculated  using  the  approximation  formulas  for  agreed  to  within  5%  with  those  obtained  by 
numerical  integration  of  Eq.  (16)  using  the  exact  Mie-theory  values  of  0,,  for  wavelengths  between  1 
/xm  and  1 1  ;txm  and  for  3.6  <  p  <  4.4,  1.2  <  n  <  1.65,  and  0  ^  k  <  0.15. 

For  calculating  extinction  coefficients  from  the  measured  aero.sol  size  distributions  and  for  investi¬ 
gating  the  contribution  of  selected  particle  size  ranges  to  the  value  of  k,.,  a  program  was  written  for  the 
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HP-9845  desktop  computer.  This  program,  which  uses  the  approximation  formulas  for  0^,  is  a  revised 
version  of  one  developed  by  Gathman  (1982)  and  was  written  to  accept  data  in  a  more  general  form. 

The  particle-si/e-distribution  data  are  put  in  analytic  form  for  use  in  the  program  by  fitting  the 
data  points  to  the  sum  of  two  or  three  log-normal  curves  of  the  form 

dN/dr  =  (17) 

by  a  computer  program  which  first  fits  a  log-normal  curve  to  the  three  (or  more)  data  points  giving  the 
counts  from  the  channels  for  the  largest  particles.  A  second  log-normal  curve  is  fitted  to  the  data 
points  remaining  after  the  first  log-normal  curve  is  subtracted  from  the  original  data  points.  If  neces¬ 
sary,  a  third  log-normal  curve  can  be  fitted  to  obtain  an  adequate  fit  to  the  smallest-particle  channels. 
The  method  used  is  that  given  by  Gathman  (1983).  The  computer  then  plots  the  sum  of  the  computed 
log-normal  curves  along  with  the  original  data  points. 

Royco-counter  data  points  do  not  always  merge  smoothly  with  mobility-analyzer  data  points, 
requiring  a  decision  as  to  how  to  weight  the  data  points  in  the  overlap  area.  Also,  when  it  became  evi¬ 
dent  that  aerosol  particles  less  than  0.1  /jim  in  radius  make  little  contribution  to  the  extinction  at 
wavelengths  greater  than  0.3  Moa,  the  best  log-normal  curves  were  selected  on  the  basis  of  the  fit  to 
points  for  a  radius  of  0.089  ^cm  or  larger.  For  some  distributions,  two  log-normal  curves  give  a  good  fit 
to  all  points  for  a  radius  greater  than  0.01  /xm,  and  occasionally  one  log-normal  curve  would  give  a  rea¬ 
sonable  fit  to  all  uata  points  (a  reasonable  fit  being  defined  as  a  fit  that  would  result  in  a  computed 
extinction  within  20%  of  the  value  computed  for  the  best  fit). 

The  largest  particle  size  used  in  obtaining  the  log-normal  fits  to  the  measured  data  points  was  a 
radius  of  2.2  /urn.  In  computing  the  extinction  from  Eq.  (16)  using  the  fitted  size  distribution,  the 
integration  was  cut  off  at  a  radius  of  3  /um.  Thus  the  computed  extinction  at  all  wavelengths  is  that  due 
to  particles  smaller  than  3  /am  in  radius.  Although  particles  larger  than  3  fxm  normally  make  a  negligi¬ 
ble  contribution  to  the  extinction  at  0.55  /am  in  wavelength,  they  can  make  an  appreciable  contribution 
to  the  extinction  at  wavelengths  greater  than  2.0  /am. 

The  real  and  imaginary  parts  of  the  index  of  refraction  of  water  are  well  known.  The  complex 
refractive  index  of  marine  and  continental  particles  is  not  as  well  known  and  is  usually  computed  on  tne 
basis  of  models  which  assume  that  atmospheric  aerosol  particles  consist  of  various  percentages  of 
water-soluble  material,  dustlike  material  (a  portion  representative  of  midlatitude  soil  conditions),  and 
soot  and  that  the  various  components  are  uniformly  mixed  with  each  other.  The  effect  of  relative 
humidity  (effect  of  water  condensed  on  the  particle)  on  the  refractive  index  is  taken  into  account  by  the 
relation  given  by  the  equations 

n  —  +  (hq  —  n,.)  (ro/r)^  (18a) 

and 

K  =  K»,  -(-  (kq  —  Kh)  (ro/r)^,  (18b) 

where  and  k„  are  the  real  and  imaginary  parts  of  the  refractive  index  of  water,  and  k,i  are  the  real 
and  imaginary  parts  of  the  refractive  index  of  the  particle  in  the  dry  state,  and  r^/r  is  the  ratio  of  the 
size  of  a  particle  in  the  dry  state  (at  0%  relative  humidity)  to  the  size  at  a  particular  humidity. 

Index-of-refraction  data  for  various  dry  aerosol  components  are  taken  from  Paltridgc  and  Platt 
(1976)  and  Shettle  and  Fenn  (1979).  Values  are  stored  for  water  and  the  following  dry  components: 
sea  salt,  water-soluble  residue,  dustlike  material,  and  soot.  Extinction  can  be  computed  for  aerosol  par¬ 
ticles  composed  of  water  and  a  selected  dry  component,  excluding  soot  (which  cannot  be  used 
exclusively  as  the  dry  component  because  the  values  of  index  of  refraction  exceed  the  limits  for  which 
the  approximation  formulas  for  have  been  shown  to  be  valid),  or  for  water  and  a  mixture  of  dry 
materials  including  soot. 
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Seven  relationships  for  the  size  of  a  particle  as  a  function  of  relative  humidity  are  available  for  use 
in  computing  the  index  of  refraction  from  Eqs.  (18).  Four  relationships  are  from  Fitzgerald  (1978): 

(r/r„)’  =  1  +  .u/d  -  W/lOO),  (19) 

where  H  is  the  relative  humidity  in  percent  and  ,e  =  1  (maritime  aerosol),  0.6  (modified  maritime),  0.3 
(modified  eonlinenuil) ,  or  0.1  (eoniinental).  1  hree  relationships  are  from  Paltridge  and  Platt  (1976): 

(/-//■„)'  =  1/(1  -  ///lOO)’  (20) 

w'here  //  is  the  relative  humidity  and  d  =  5.8  (continental),  4.4  (normal  or  mixed),  or  3.9  (maritime). 

Extinction  is  computed  at  selected  wiavelengths,  with  the  program  being  especially  set  up  for  com¬ 
putations  at  0.55  /cm.  The  extinction  at  selected  wavelengths  is  stored  in  memory  during  computation 
and  plotting  of  the  extinction  and  can  be  printed  after  plotting  is  finished.  Extinction  due  to  a  selected 
particle-size  range  can  also  be  computed  for  selected  wavelengths  as  well  us  being  computed  and  plotted 
for  all  wavelengths. 

Since  the  computed  extinction  may  vary  by  25'’/i  or  more  between  different  log-normal  fits  to  the 
same  data,  depending  on  the  weighting  of  points  in  the  overlap  region  of  the  Royco-counter  and 
mobility-analyzer  data,  variations  in  extinction  of  several  percent  due  to  uncertainties  in  the  index  of 
refraction  of  the  particles  are  not  significant.  Thus  the  computed  value  of  extinction  is  relatively  insen¬ 
sitive  to  whether  the  growth  relationships  of  Fitzgerald  or  Paltridge  and  Platt  are  used  in  Egs.  (18)  or  to 
the  choice  of  aerosol  type  (maritime  versus  modified  maritime,  etc.).  For  relative  humidities  less  than 
about  80'’/o,  the  effect  of  a  lO'/n  error  in  the  measured  relative  humidity  on  the  value  of  the  index  of 
refraction  results  in  only  a  few  percent  change  in  computed  extinction. 

Computation  of  the  extinction  coefficient  for  curve  1  in  the  left-hand  plot  of  Fig.  41  illustrates  the 
method.  Figures  8  and  9  show  two  different  log-normal  fils  corresponding  to  curve  I  of  Fig.  41.  Fig¬ 
ure  8  represents  a  good  log-normal  fit  to  all  data  points  of  radius  greater  than  0.5  /im,  and  Fig.  9  shows 
a  best  fit  when  the  smallest-particle  Royco  channel  is  lightly  weighted.  (The  agreement  between  meas¬ 
ured  and  computed  extinction  is  usually  improved  if  the  smallest-particle  Royco  channel  is  lightly 
weighted.  The  fit  to  the  remaining  points  is  also  improved.)  Figures  10  and  11  are  plots  of  extinction 
versus  wavelength  for  the  log-normal  fits  shown  in  Figs.  8  and  9  respectively.  The  upper  curve  in  each 
figure  is  the  total  extinction,  and  the  other  curves  are  extinctions  due  to  particles  m  selected  size 
ranges. 

The  upper  and  lower  sections  of  Table  1  show  the  computed  values  of  the  extinction  cocffiLient 
for  the  log-normal  fits  in  Figs.  8  and  9  respectively.  Tabulated  is  the  total  extinction  coefflccnt  (bot¬ 
tom  line  in  each  section)  and  the  partial  extinction  due  to  particles  in  selected  size  ranges  Fi  this  case 
the  different  log-normal  fits  to  the  same  data  give  rise  to  a  28%  difference  in  the  total  coefficient  The 
largest  differences  in  partial  extinction  are,  as  would  be  expected,  in  the  size  range  0  2  to  0  4  /im, 
which  includes  the  smallest-particle  Royco  channel  at  0.24  //m.  The  extinction  by  particles  below  0  1 
/im  is  very  small,  although  these  particles  are  most  of  the  particles  counted 


Radon  Instrumentation 

Emanations  of  the  radioactive  rare  gas  radon  (”^Rn)  from  the  earth's  surface  and  subsequent 
dispersal  in  the  atmosphere  are  of  interest  in  the  study  of  vertical  diffusion  (Bradley  and  Pearson,  1970, 
and  Hosier,  1969),  in  the  study  of  dust  transport  (Prospero  and  Carlson,  1970),  and  as  an  indicator  of 
whether  or  not  an  air  mass  had  a  trajectory  recently  over  or  near  a  land  mass 
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Fig.  8  —  An  example  of  the  best  fit  of  three  log¬ 
normal  functions  to  all  experimental  points  on  the  size 
distribution  shown  as  curve  I  in  Fig  40,  The  crosses 
are  experimental  data,  and  the  dashed  line  is  the  sum 
of  the  three  log-normal  curves 


RADIUS  i/im) 

Fig.  9  —  Same  example  as  Fig  8  except  that  the  smal¬ 
lest  channel  of  the  optical  counter  was  ignored  in  cal¬ 
culating  the  best  fit 


Fig  10  —  Fxiinction  as  a  function  of  wavelength  for  the 
fitted  size  distribution  shown  in  Fig  8 
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Table  1  —  Computed  values  of  the  extinction  coefficient 
for  the  log-normal  curve  fits  shown  in  Fig.  8  and  9.  The 
measured  relative  humidity  was  14%.  The  panicle  index 
of  refraction  at  74"/i  relative  humidity  was  computed  from 
Eqs.  (18)  and  (19)  using  ^  =  1. 


Ra()iu5s 

Range 

(fim) 

Exiinction  Coeffieieni  (10''*/m) 
at  Different  Wavolenglh.s  (A» 

0  4 

0  3  fi  m 

0  33  litn 

0.63  fi  m 

1  /uni 

2  ^ni 

Computed  for  the  Fit  in  Fig  8  and  Plolled  in  Fig  10 

0.04-0.07 

0016 

OOOS 

0007 

0007 

000 

0  00 

0.07-0  10 

0.068 

0,028 

0023 

0  021 

001 

0  00 

0.10-0  20 

0,6.VS 

0  333 

0  271 

0  189 

0  03 

0  01 

0.20-0.40 

0  671 

0,504 

0.427 

0.3.30 

oil 

0  02 

0.40-0.60 

0  16.3 

0  214 

0.223 

0224 

0  1.3 

0  0.3 

0.60-0.80 

0  107 

0  123 

0  1.39 

0  171 

0  18 

0(13 

0.80-1.00 

0.078 

0.079 

0  079 

0  088 

0  14 

006 

1.00-3.00 

0,128 

0  128 

0  123 

0  123 

0.17 

017 

0.04-3.00 

1.808 

1  387 

L 

1  095 

0  74 

0.31 

Computed  for  the  Fit  in  Fig  9  and  Plotted  in  Fig  1 1 

0.04-0  07 

0016 

0.008 

0007 

0  007 

OOO 

0  00 

0.07-0.10 

0068 

0.023 

0023 

0  021 

001 

000 

0  10-0  20 

0.743 

0,402 

0.307 

0  213 

006 

001 

0.20-0.40 

1  104 

0  841 

0  716 

0.337 

0  19 

0.0.3 

0  40-0.60 

0.213 

0.271 

0  2.30 

0.276 

0  13 

0  03 

060-0,80 

0  100 

0  1 1 3 

0  1.30 

0  160 

0  16 

004 

0  80-1  00 

0.07f 

0  073 

0073 

0(184 

0  14 

0  06 

1,00-3  no 

0  123 

0.124 

0  1 24 

0  124 

0  16 

0  17 

0  04-3  00 

2  341 

1,797 

1  607 

1  ,367 

0  82 

0  32 
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Radon  is  part  of  the  uranium  series,  which  starts  with  and  ends  with  stable  ’"'’Pb.  '-’Rn  is 
the  daughter  of  ’’’’Ra.  The  source  of  atmospheric  radon  is  the  trace  amounts  of  radium  present  m  all 
soils  in  varying  concentrations.  Gaseous  --’Rn  diffuses  from  the  earth's  surface  essential!)  free  of  an> 
daughter  activity  and  reaches  secular  equilibrium  with  its  daughter  products  through  ’''’Po  within  hours 

Radon  daughter  products  are  electrically  charged  and  chemically  reactive.  These  particles  rapidly 
become  attached  to  natural  aerosol  particles  and  can  be  collected  on  filters.  Short-lived  radioactive 
atoms  can  live  out  their  lives  in  a  free  condition  in  clear  air,  but  it  has  been  demonstrated  (Saunders. 
Patterson,  and  Lockhart,  1968)  that  these  unattached  radon  daughters  are  collected  by  cellulose-glass- 
fiber  and  glass-fiber  filters  with  nearly  the  same  efficiency  as  aerosol  particles  are  collected.  The  contri¬ 
bution  of  long-lived  ^'®Pb  and  subsequent  daughter  products  in  the  uranium  series  to  activity  collected 
on  the  filter  has  been  shown  to  be  negligible  (Blifford,  Lockhart,  and  Rosenstock,  1952). 

A  typical  sample  collector  consists  of  a  filter  holder  and  a  blower  to  draw  air  through  the  filter.  A 
detector  consists  of  a  plastic  scintillator  optically  coupled  by  silicone  grease  to  a  photomultiplier  tube. 
The  plastic  scintillator  is  covered  with  0.05  mm  of  aluminum  to  absorb  alpha  particles  and  provide  a 
light  shield.  Total  beta  activity  on  the  filter  changes  slowly  with  age  after  about  an  hour.  Plastic  scintil¬ 
lators  0.25  mm  thick  were  selected,  because  they  have  been  shown  to  be  thick  enough  to  give  adequate 
light  output  from  passage  of  beta  rays  (Harley  et  al.,  1962)  and  yet  thin  enough  to  be  essentially  tran¬ 
sparent  to  gamma  rays  and  have  a  low  background  counting  rate. 

The  first  radon  counter  constructed  in  the  NRL  Atmospheric  Physics  Branch  was  constructed  to 
monitor  the  maritime  air  at  a  Barbados  field  site  for  possible  contamination  by  continental  air  that  is 
frequently  transported  across  the  Atlantic  from  Africa.  Commercial  instrumentation  did  not  have  the 
required  sensitivity  and  portability.  The  first  instrument  combined  the  radon  measurement  techniques 
developed  by  the  Chemistry  Division  (Lockhart  et  al.,  1965)  with  sensitive  beta-particle  detectors 
developed  by  the  Radioactivity  Branch  of  the  Nucleonics  Division  in  response  to  a  need  far  monitoring 
the  atmosphere  of  the  USS  Nautilus  for  fission  products  (anonymous,  1960). 

An  improved  version  (Larson,  1973)  was  used  by  personnel  of  the  Atmospheric  Physics  Branch 
and  the  Chemical  Oceanography  Branch  on  aircraft  and  ships  to  measure  radon  in  remote  areas  and  to 
relate  the  radon  to  concentrations  of  other  materials  of  continental  origin  and  to  atmospheric  processes. 
This  instrumentation  could  measure  radon  down  to  a  few  tenths  of  a  picocurie  per  cubic  meter  and  thus 
was  about  two  orders  of  magnitude  more  sensitive  than  the  available  commercial  instruments,  which 
were  designed  for  monitoring  continental  radon  levels  of  about  100  pCi/m\  or  the  radioactivity  in 
homes  or  uranium  mines,  where  the  radon  is  much  higher. 

The  success  of  this  instrumentation  prompted  the  building  of  the  Automatic  Radon  Counter 
(ARC)  (Larson  and  Bressan,  1978)  (U.S.  Patent  4,140,912)  and  enabled  automatic  and  frequent  collec¬ 
tion  of  data.  Thus  data  can  be  collected  around  the  clock,  even  when  it  is  unsafe  for  personnel  to  go 
out  on  deck  as  required  by  the  nonautomated  system. 

The  NRL  Environmental  Chemistry  Branch  has  combined  the  automatic  radon  measurement  with 
an  airborne  salt  measurement  using  the  same  filter  papers.  This  instrumentation  (Bressan.  1981 1  has 
enabled  unattended  round-trip  shipboard  measuring  of  radon  and  salt. 

The  Atmospheric  Physics  Branch  has  recently  used  a  simplified  version  called  the  .Mmost 
Automatic  Radon  Counter  (AARC).  The  ARC  had  a  tape  drive  which  increment-  the  filter  paper  from 
collection  position  to  counting  position  after  each  sample  is  collected  and  thus  has  a  cl '  m  filter  for  each 
sample.  Samples  can  be  collected  as  often  as  every  5  or  10  minutes  (with  sonve  loss  of  sensitivity  I. 
The  AARC  uses  the  same  filter  for  successive  samples.  The  blower  and  counter  ,ire  automatic.illy 
cycled  to  enable  relatively  unattended  operation,  but  the  filter  must  be  changed  every  1  to  .1  days,  and 
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Ihc  t'ri.'i|ucnLy  ol'  samples  can  be  no  more  than  one  per  hour.  The  ad\aniage  of  the  AARC  over  the 
AR(  IS  greatly  reduced  cost  and  size  of  those  components  that  must  be  exposed  to  the  elements. 


The  radon  data  presented  in  this  report  were  taken  with  the  .AARC  version  and  are  used  primarily 
as  a  diagnostic  tool  to  help  identify  the  source  of  air  masses. 

Supporting  Instrumentation 

Atmospheric  research  requires  measurements  in  undisturbed  air  free  of  local  sources  of  pollution. 
Also  required  is  shelter  for  personnel  and  instrumentation,  and  air  conditioning  may  be  required  for 
sophisticated  instrumentation  and  computers.  As  the  number  and  complexity  of  our  instruments 
increased,  it  became  increasingly  difficult  to  find  existing  locations  with  adequate  exposure  for  air  sam¬ 
pling  and  vents  and  preferably  outside  the  ship's  aerodynamic  boundary  layer.  The  time  required  to  set 
up  and  dismantle  instrumentation  on  site  or  aboard  ships  became  increasingly  costly. 

To  alleviate  these  recurring  problems,  it  was  necessary  to  purchase  a  movable  equipment  shelter 
which  could  be  transported  to  the  desired  field  sites.  The  shelter  had  to  be  large  enough  for  all  the 
present  and  immediately  foreseeable  instrumentation  and  operating  personnel,  had  to  be  portable  for 
easy  transportation  to  field  sites  and  for  relocation  at  sites,  and  had  to  be  small  enough  for  installation 
on  the  flying  bridge  of  the  USNS  Hayes  or  other  research  ships. 

The  selected  design  was  a  mobile  equipment  shelter  2.5  m  wide  and  high  and  5.5  m  long 
manufactured  by  Ekto  Corporation  of  Sanford,  Maine,  This  unit  was  selected  because  it  is  strong 
enough  to  withstand  frequent  deployment,  is  light  enough  to  be  loaded  aboard  ship,  and  was  competi¬ 
tive  in  cost.  The  shelter  was  acquired  with  a  triple-axle  transporter  for  trailing  it  to  a  field  site.  Corner 
jacks  enable  stabilizing  and  leveling  at  field  sites.  The  shelter  can  easily  be  removed  from  the  tran¬ 
sporter  for  shipboard  use. 

The  200-ampere  service  in  the  shelter  is  more  than  sufficient  for  air  conditioning  and  existing 
instrumentation.  There  are  five  data-cable/air-intake  ports  and  an  extra-wide  door  on  one  end  for  load¬ 
ing  equipment.  The  top  is  flat,  has  a  removable  safety  railing,  and  is  strong  enough  to  support  a  tower 
to  further  elevate  instrumentation.  A  large  window  in  one  side  and  small  windows  in  each  door  provide 
visibility  of  external  conditions.  The  shelter  has  built-in  lift  and  tie-down  points,  adequate  to  secure  it 
during  rough  seas. 

The  shelter  has  a  small  area  for  a  desk  and  has  rails  attached  to  the  floor  for  securing  instrument 
racks.  Most  of  our  equipment  can  be  installed  and  checked  out  and  altered  if  necessary  at  NRL.  thus 
greatly  reducing  setup  time  on  site.  Once  equipment  has  been  installed  in  the  shelter  for  one  field 
experiment,  reinstallation  time  for  subsequent  field  work  is  minimal. 

In  addition  to  the  instruments  previously  discussed,  the  shelter  housed  meteorological  instrumen¬ 
tation  for  the  measurement  of  temperature,  dewpoint,  wind  speed,  and  wind  direction.  The  sensors  for 
these  measurements  were  mounted  on  the  roof  of  the  shelter. 

An  instrument  which  we  refer  to  as  the  variable-humidity  nephelometer  has  been  developed  at 
NRL  and  was  successfully  used  during  the  Washington,  D.C.,  deployment.  It  was  also  taken  on  the 
//ayes  cruise.  This  instrument  in  its  present  design  is  manually  operated  and  requires  full  attention  of 
at  least  one  person  during  its  operation.  It  was  not  operated  during  the  cruise  becau.se  of  manpower 
shortage  and  enhanced  awkwardness  of  operation  due  to  ship's  motion.  This  instrument  and  data 
obtained  during  the  Washington,  D.C.,  deployment,  are  described  by  Fitzgerald  et  al.  (1982). 
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RESULTS 


In  early  June  W?*)  the  instrument  shelter  was  plaectl  near  the  NRL  waterfront,  and  data  were 
taken  for  a  number  of  days  to  test  equipment  and  develop  procedures.  From  there  the  equipment  was 
transported  to  Wallops  Island,  Virginia,  and  placed  on  one  of  the  rocket  launching  pads  about  80  m 
from  the  coast  (at  high  tide).  Data  were  taken  at  this  site  from  18  June  until  29  June  1979.  The 
equipment  was  returned  to  NRL,  and  additional  data  were  collected  during  July  1979.  In  the  fall  of 
1979  the  opportunity  to  participate  with  the  NRL  Marine  Biology  Branch  on  a  cruise  scheduled  early  in 
the  following  year  presented  itself.  The  USNS  Hayes  finally  cruised  in  April  1980,  sailing  from 
Norfolk,  Virginia. 

The  measurements  presented  in  this  report  are  those  taken  on  these  three  deployments.  Many 
data  were  obtained,  and  no  attempt  will  be  made  to  present  all  of  the  data  and  required  analysis. 
Instead  selected  periods  will  be  chosen  for  presentation.  The  greatest  emphasis  will  be  on  the  shipboard 
measurements.  This  emphasis  reflects  the  particular  interest  of  the  Navy  in  marine  aerosols  and  the 
effect  they  have  on  naval  systems. 


Washington,  D.C.,  Measurements 

Although  the  aerosol  program  was  initiated  to  look  at  the  marine  and  coastal  environment,  the 
development  and  testing  of  the  instrumentation  at  NRL  gave  us  an  opportunity  to  look  at  urban  aero¬ 
sols  as  well.  Some  samples  of  these  data  will  be  presented  to  contrast  the  urban  environment  with  the 
marine  environment.  Also,  an  air  mass  advecting  from  the  U.S.  east  coast  over  the  Atlantic  may  have 
passed  over  some  urban  areas  during  its  transit,  and  the  urban  aerosol  will  contribute  to  the  history  of 
the  aerosol  load  in  that  air  mass.  Furthermore  study  of  the  urban  aerosol  will  contribute  to  our  under¬ 
standing  of  the  physical  processes  affecting  the  evolution  of  the  size  distribution. 

The  size  distribution,  scattering  coefficient,  and  meteorological  variables  were  measured  routine's 
during  the  day  between  19  July  and  30  July.  On  several  occasions  measurements  were  also  made  c.o 
ing  the  night.  Radon  data  were  not  taken  at  the  NRL  site,  because  its  greatest  utility  is  for  distingqmi- 
ing  marine  air  masses  from  those  which  have  recently  been  over  land  in  a  coastal  and  marine  environ¬ 
ment.  Data  from  the  variable-humidity  nephelometer  were  obtained  only  at  the  NRL  site.  The  hygros- 
copicity  parameter  S"  of  small  particles  was  measured  on  a  number  of  occasions  using  the  method 
described  earlier.  The  results  of  these  hygroscopicity  measurements  will  be  presented  later  in  this 
report  along  with  results  from  other  locations. 

As  is  well  known,  the  urban  environment  contains  a  large  number  of  small  particles  generated  by 
man’s  activity.  Figure  1  shows  a  rather  typical  urban  size  distribution  (average  of  four  individual  distri¬ 
butions)  taken  midmorning  on  July  30  and  labeled  Washington,  D.C.  During  this  particular  time  the 
scattering  coefficient  bj  was  about  2.3  x  10  Vm,  the  wind  was  from  the  north-northeast,  and  the  total 
number  of  particles  was  about  33,000/cm’.  Total  counts  at  the  site  during  the  2-week  interval  were 
typically  in  the  range  10,000/cm’  to  40,000/cm’. 

Figure  12  shows  size  and  area  distributions  during  two  different  meteorological  conditions  with 
very  different  visibilities.  The  size  and  area  distributions  labeled  1  are  each  the  average  of  five  distribu¬ 
tions  during  an  interval  on  23  July  when  conditions  remained  fairly  stable  with  good  visibility 
(bj  =  0,7  X  10  '’/m)  and  a  total  aerosol  count  of  15,300  /cm’.  The  wind  was  from  the  south,  the  sky 
was  mostly  cloudy,  and  the  relative  humidity  was  about  80%.  A  high-pressure  area  off  the  mideast 
Atlantic  coast  (Bermuda  high)  was  causing  air  to  flow  irom  the  Atlantic,  over  a  portion  of  North  Caro¬ 
lina  and  Virginia,  and  into  the  Washington,  D.C.,  area. 
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The  size  and  area  distributions  labeled  2  are  each  the  average  of  five  distributions  taken  on  the 
afternoon  of  30  July.  The  visibility  was  poor  (b^  =  3.6  x  10  Vm),  the  total  particle  count  was  about 
23,000/cm^,  and  the  relative  humidity  was  about  77%.  The  wind  was  light  from  the  north-northeast 
and  resulted  from  a  northerly  flow  around  a  weak  high-pressure  area  centered  over  northwest  Pennsyl¬ 
vania. 


Area  distributions  similar  (at  radii  greater  than  0.1  /xm)  to  the  one  measured  on  30  July  were  also 
obtained  on  20  July  (not  shown  in  the  plots).  On  20  July  Washington  was  under  a  stagnant  high- 
pressure  ridge  extending  from  Illinois  to  New  York  and  centered  over  Pennsylvania,  visibilities  were 
also  low  (6,  =  3.2  x  10  Vm),  and  the  relative  humidity  was  65%. 

Since  the  measured  scattering  coefficients  associated  with  curves  1  and  2  of  fig.  12  are  very 
different,  it  is  of  interest  to  compare  the  calculated  and  measuicd  values  for  these  times.  Figures  13 
and  14  show  the  extinction  coefficients,  calculated  from  the  size  distributions  by  the  method  de,'>cribed 
earlier,  for  curves  1  ana  2  respectively.  The  top  curve  of  Figs.  13  and  14  is  the  total  extinction.  The 
other  curves  give  the  contributions  of  various  size  classes.  In  both  cases  particles  in  the  radius  range 
0.1  to  0.3  /am  are  the  largest  contributors  to  extinction  at  short  wavelengths.  However,  between  0.5 
/xm  and  2.0  /am  in  wavelength  the  relative  contribution  of  the  different  size  classes  in  the  two  cases 
varies  significantly.  Table  2  compares  the  measured  .scattering  coefficient  and  the  extinction  coefficient 
calculated  from  the  size  distributions. 


Table  2  —  Comparison  of  measured  and 
calculated  extinction  coefficient 


Extinction  (10  Vm) 

Date 

Measured 

Calculated 

0.55  um 

0.55  ,xm 

0.50  ,xm 

0.63  ,xm 

23  July 

0.70 

0.80 

0.89 

0.68 

30  July 

3.6 

4.50 

4.67 

4.39 

23 
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Fig.  13  —  FxtinOiitn  cakul.ilions  basecl  on  the  si^c  tlistribulion  shown 
in  curve  I  in  the  left-hand  plot  of  Fig  13 
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Fig  14  —  I  xtinOion  c.ilcul.itions  hasevi  on  the  si/e  disirihulion  shown 
in  curve  3  in  the  left-hand  plot  of  Fig  I  3 
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Wallops  Island  Measurements 

The  Site 

After  completion  of  a  test  period  at  NRL,  the  equipment  was  moved  to  NASA's  Wallops  Island 
rocket-launching  site.  The  instrument  shelter  was  placed  on  one  of  the  concrete  launching  pads  80  m 
from  a  sand  beach.  The  site  itself  is  about  5  nmi  from  the  main  NAS.A  base  and  tracking  center  (Fig. 
15)  and  is  surrounded  by  marshland.  To  the  east  is  the  Atlantic  Ocean,  immediately  to  the  west  is  the 
Delmarva  Peninsula,  and  farther  to  the  west  is  the  Chesapeake  Bay.  The  Delmar\a  Peninsula  is  rela¬ 
tively  flat  farmland  with  no  significant  local  sources  of  industrial  pollution. 


'  I'H:.  I-''  —  W.illdps  lsl.(ntl  rnc.isiiromcnl  ^ll^;  I  ho  mo.iMiromonl  Ido.ilidi)  m.iikoi)  hi  x 

[ 

j  The  coastal  site  was  selected  to  study  the  difference  in  aerosols  in  air  masses  of  continental  origin 

I  and  of  marine  origin.  Air  masses  arriving  at  the  site  from  the  west  should  represent  the  rural  east-coast 

I  environment,  whereas  air  masses  arriving  from  the  east  would  presumably  represent  a  more  marine 

environment.  Air  arriving  from  the  east  would  pass  over  the  coastal  surf  zone  and  possibly  be  modified 
5  by  surf-generated  particles.  In  addition  a  range  of  conditions  between  pure  continental  and  pure  mari- 

i  time  would  result  from  air  masses  whose  recent  trajectories  had  passed  over  both  marine  and  continen- 

I  tal  regions. 

1 

The  sample  air  for  the  mobility  analyzer  and  an  Environment/One  condensation-nuclei  counter 
was  drawn  from  about  4.5  m  above  ground  through  a  25-mm-l.D.  tube  at  a  volume  flow  rate  of  about 
0.07  m^/min.  Air  for  the  nephelometer,  also  housed  in  the  shelter,  had  its  own  38-mm-l.D.  sampling 
duct.  The  sensor  of  the  Royco  optical  particle  counter  was  mounted  on  top  of  the  shelter  together  with 
the  radon  monitor  and  the  meteorological  instrumentation. 
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Meteorological  Conditions 

The  synoptic  conditions  over  the  eastern  United  Stales  during  the  deployment  were  sufficiently 
varied  that  we  were  able  to  obtain  aerosol  samples  in  air  masses  with  a  variety  of  histories  Hstimates 
of  the  recent  histories  of  the  air  being  measured  were  based  on  NOAA  surface  weather  charts,  local  and 
area  wind  patterns,  and  radon  concentrations.  Local  meteorological  conditions,  weather  forecasts  and 
radiosonde  data  were  provided  by  the  NOAA  weather-service  support  facility  at  Wallops  Flight  Center. 

Radon  concentrations,  CCN  measurements  at  0.7%  supersaturation,  and  local  winds  at  Wallops 
Island  are  shown  in  Fig.  16.  (The  times  indicated  by  arrows  along  the  top  border  are  for  use  in  a  later 
discussion.)  The  radon  concentrations  of  between  2  pCi/m^  and  3  pCi/m’  on  27  and  28  June  are  low 
and  indicate  that  the  air  had  a  long  trajectory  over  water.  We  will  refer  to  this  type  of  air  mass  as  being 
marine.  The  higher  concentrations  near  100  pCi/m^  indicate  a  recent  history  of  continental  air. 


JUNE  1979 
ibl 
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CCN  concentrations  are  plotted  along  with  radon  because  it  has  been  observed  that  radon  and 
CCN  correlate  well  at  coastal  and  oceanic  locations  (Larson  et  al.,  1979).  The  source  strength  of  an> 
land  area  will  be  different  for  radon  and  CCN,  since  the  radon  source  strength  depends  on  the  radium 
content  and  physical  nature  of  the  soils  and  the  CCN  source  strength  is  related  to  natural  and  anthropo¬ 
genic  gas  and  particle  sources. 

The  first  samples  were  taken  late  on  18  June  in  continental  air  behind  a  cold  front  that  had  moved 
through  the  area  earlier  that  day.  By  the  morning  of  19  June,  '.he  radon  concentration  had  decreased, 
and  the  winds  and  synoptic  meteorology,  combined  with  the  lower  radon  and  CCN  readings,  suggested 
some  subsidence  of  cleaner,  high-altitude  air  within  the  high-pressure  area,  centered  over  the  Great 
Lakes.  By  late  afternoon,  the  winds  were  from  the  northeast,  and  some  maritime  inllucnce  had  mixed 
into  the  air. 

The  synoptic  situation  on  20  June  was  such  that  the  air  movements  were  not  well  defined.  .A 
high-pressure  area  over  New  York  with  a  weak  pressure  gradient  encouraged  a  noriherlt  flow  ol  more 
continental  air  into  the  area,  but  a  weak  low-pressure  perturbation  over  North  Carolina  brought  mari¬ 
time  air  into  the  area  by  afternoon.  On  21  June  a  weak  high-pressure  ridge  over  the  east  coast  of  the 
United  States  brought  a  mixture  of  maritime  and  continental  air  into  the  area  and  possibly  even  some 
industrial  contamination  from  the  Norfolk,  Virginia,  area. 

The  air  became  mostly  continental  on  22  June  and  remained  that  wav  through  24  June.  Winds 
were  generally  southwesterly  on  22  and  23  June  and  northerly  and  northeasterly  on  24  June  Some  rain 
occurred  during  this  period. 

Throughout  the  3-day  interval  of  mostly  continental  air,  the  correlation  coefficient  between  radon 
and  the  height  of  the  inversion  was  -0.8.3,  indicating  a  good  inverse  correlation  Hsu  et  al.  (1980) 
reported  a  correlation  of  -0.96  for  radon  and  inversion  height  on  the  gulf  coast  of  Florida  near  Panama 
City  in  continental  air  moving  off  the  coast  from  the  central  United  .States  .A  good  inverse  correlation 
between  radon  and  inversion  height  indicates  that  the  air  in  the  mixed  laver  had  been  of  relatively  uni¬ 
form  continental  origin  over  the  preceding  day  or  so. 

The  interval  25  through  28  June  was  charactcri/ed  by  a  transition  from  continental  air  e.irlv  on  25 
June  to  maritime  air  late  on  27  June  The  air  was  mostly  maritime  on  28  June,  but  the  maritime  flow 
was  overcome  by  local  offshore  flow  early  on  28  June,  which  introduced  a  continental  contamination  to 
the  maritime  air.  The  synoptic  transition  as  shown  bv  the  surface  charts  is  shown  in  I  ig  I  7  On  2s 
June  the  clockwise  circulation  around  a  high-pressure  area  centered  over  the  (ireat  Lakes  at  0500  LOT 
was  bringing  mostly  continental  air  to  the  site.  By  late  in  the  afternoon,  the  winds  became  onshore 
(easterly),  and  the  radon  concentration  decreased  to  the  teens,  indicating  a  mixture  of  continental  and 
marine  air. 

By  26  June,  the  high  pressure  was  centered  over  New  York  State,  and  the  continental  char.icter  of 
the  air  continued  to  decrease.  The  air  reaching  the  site  had  spent  about  2  days  over  the  ocean  and  anv 
near  surface;  thus  high-radon  air  was  uniformly  distributed  throughout  the  mixed  laver  and  mav  also 
have  been  diluted  by  subsident  air  or  by  additions  of  maritime  air  By  27  June,  the  high  pressure 
center  was  off  Cape  Cod,  but  the  air  flow  at  the  site  remained  easterly,  and  the  radon  concentrations 
indicate  maritime  air  by  radon  levels  that  arc  near  background  levels  for  the  northern  North  Atlantic 
On  28  June  much  of  the  eastern  seaboard  was  under  the  influence  of  a  large  high  pressure  with  a  low 
gradient  of  pressure. 

The  correlation  between  radon  and  CCN  was  good  the  second  6  days  (lower  plot  in  Fig  161,  as 
would  be  expected  in  air  that  had  approximately  similar  origin  but  spent  an  increasinglv  greater  time 
away  from  land.  The  correlation  was  0.82  The  correlation  coefficient  between  radon  and  CCN  was 
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Fig  1 7  —  Movcmcni  of  a  high-pressure  s.\  slem  across 
the  northern  Unitet)  .States 


0.72  for  the  entire  12  days  of  operation  at  Wallops  Island.  The  lower  overall  correlation  is  a  conse¬ 
quence  of  the  different  source  areas  for  different  days. 

Size  Distributions 

Size  distributions  were  taken  daily  throughout  the  12  days  except  on  24  June.  There  were  periods 
during  the  morning  hours  when  no  data  were  obtained.  The  mobility  analyzer  would  malfunction,  and 
much  time  would  be  spent  testing  for  nonexistent  leaks  in  the  filtered  sheath  air  system.  It  was  finally 
determined  that  the  source  of  the  problem  was  turbulent  convection  in  the  analyzer  arising  from  tran¬ 
sient  temperature  differences  during  warmup  after  a  cool  night.  The  problem  was  simply  remedied  by 
leaving  the  sheath-flow  pump  running  continuously,  maintaining  the  entire  unit  at  the  same  tempera¬ 
ture  overnight.  Even  though  the  system  was  automated,  data  were  not  taken  continuously  throughout 
the  night,  because  there  was  a  problem  with  condensation  in  the  segmented  chamber  which  required 
attention  every  few  hours.  Before  the  cruise  in  April  1980  this  problem  was  cot  ected. 

Figure  1  illustrates  differences  in  the  size  distribution  obtained  at  Wallops  island  between  periods 
when  the  air  was  definitely  from  the  continent  and  when  the  air  arrived  at  the  site  after  a  long  trajectory 
over  water.  The  1 1  smaller  channels  are  from  data  taken  with  the  mobility  analyzer,  and  the  four  larger 
channels  are  from  data  taken  with  the  Royco  optical  particle  counter.  In  the  overlapping  channel  it  was 
nearly  always  the  case  that  the  mobility  analyzer  recorded  more  particles  than  the  optical  part'cle 
counter.  The  continental  size  distribution  shown  in  Fig.  1  was  taken  on  18  June,  when  the  visibility 
was  the  lowest  recorded  (b,  =  4.0  x  10  Vm)  during  our  12  days  at  Wallops  Island. 

Since  the  time  required  to  obtain  a  single  size  distribution  is  approximately  15  minutes,  any 
change  in  the  absolute  or  relative  distribution  during  that  time  would  lead  to  a  distortion  of  the  size  dis¬ 
tribution.  To  help  evaluate  Just  how  significant  natural  variations  are  during  a  measurement  interval,  it 
was  useful  to  average  a  number  of  consecutive  data  sets  and  calculate  the  size  distribution  from  the 
average  data.  To  evaluate  the  influence  of  natural  changes,  the  size  distribution  was  also  calculated 
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I'rom  the  average  data  plus  one  standard  deviation  and  (rom  the  average  data  minus  one  standard  devia¬ 
tion.  Figure  18  shows  the  result  tor  the  continental  si/e  distribution  shown  in  Fig.  1.  The  average  dis¬ 
tribution  is  calculated  from  five  data  sets  taken  over  about  2.5  hours.  The  width  of  the  data  bars 
represents  the  width  of  the  size  channel.  The  uppermost  point  on  the  vertical  bar  is  the  size  distribu¬ 
tion  calculated  from  the  average  data  plus  a  standard  deviation,  whereas  the  lowest  point  is  that  calcu¬ 
lated  from  the  average  minus  a  standard  dev  lation.  The  vertical  bars  are  not  error  bars  associated  with 
the  instrument  but  rather  are  an  indication  of  atmospheric  variability  during  that  particular  measure¬ 
ment  period.  The  variations  shown  in  Fig.  18  are  typical  in  that  the  variations  observed  in  the  smaller- 
particle  channels  are  nearly  always  greater  than  in  the  larger-particle  channels. 


001  '  . 

0.001  0.01  0.1  1  2  4  6  10 

RADIUS  Ifiml 

I  1^  IK  —  \\ciJpk’  si/c  ilistnbiitiMn  troni  live  U.il.i  sets  t.iken  at 
the  W.illnps  Kl.inJ  site  during!  a  J inieiv.il  f  he  \eitKa! 
b.irs  sh<*w  the  st.iml.nd  ile\i.ilH>n  in  the  si/c  distnhulmn  due  ti’ 
fe.d  v.iff.ftKMis  in  the  si/e  <l»slf ibiiiii>n  dti/nu’  the  inien.il 


29 


IIOI’I’I  I  MI/C.IRAII).  \M)1  \kS()N 


In  most  of  the  size  distributions  there  is  an  inflection  (change  of  slope)  in  the  optical-counter  data 
in  the  size  range  0.3  to  1.4  ^m.  In  this  size  range  there  is  a  double-valued  region  in  the  calibration 
curve  of  the  Royco  optical  particle  counter.  This  double-valued  region  is  very  sensitive  to  the  refrac¬ 
tive  index  of  the  particle.  An  effort  was  made  to  adjust  the  channel  settings  such  that  the  double¬ 
valued  region  was  entirely  within  a  single  channel.  However,  we  cannot  be  certain  that  this  was  the 
case,  since  the  refractive  index  of  natural  aerosols  is  not  well  known.  The  inflection  in  the  Royco  size 
distribution  may  result  in  part  from  particles  which  belong  in  one  size  channel  being  assigned  to  the 
other  channel. 

In  Fig.  1  the  marine  size  distribution  taken  at  Wallops  Island  indicates  a  relatively  large  number  of 
particles  in  the  third  and  fourth  Royco  channels.  Whenever  the  wind  was  from  the  ocean,  accompanied 
by  significant  surf  activity,  there  was  a  large  variability  and  increase  in  the  number  of  particles  in  the 
radius  range  above  0.5  /zm  due  to  surf-generated  particles.  Hence  it  is  doubtful  that  data  taken  at  any 
coastal  site  can  truly  represent  the  open-ocean  environment. 

During  the  12-day  period  at  Wallops  Island,  over  150  valid  size  distributions  were  obtained. 
Interpreting  the  hour-to-hour  changes  in  the  distribution  is  difficult.  Both  the  long-range  transport  and 
local  conditions  influence  the  changes.  Even  though  the  coastal  situation  offers  a  good  opportunity  to 
contrast  continental  and  marine  conditions  in  some  well-defined  meteorological  cases  (Fig  1).  at  other 
times  the  coastal  environment  is  very  complex.  Often  it  is  difficult  to  determine  whether  or  not  the 
recent  trajectory  has  been  over  land  or  water  and  to  determine  the  influence  of  local  and  sea  circula¬ 
tions  and  the  source  strength  of  surf-generated  particles.  Therefore  only  a  selective  sanii'ling  of  data  is 
presented  here. 

Several  attempts  were  made  to  observe  an  increase  in  the  number  of  small  particles  associated 
with  an  increased  production  due  to  photochemical  activity  at  sunrise,  but  no  clear  increase  was  ever 
observed.  Also,  during  the  measurements,  diurnal  effects  were  much  weaker  than  variations  caused  by 
changes  in  synoptic  conditions  and  therefore  were  not  observed. 

The  most  interesting  data  interval  occurred  between  25  June  and  28  June.  On  25  June  the  clock¬ 
wise  circulation  around  a  high-pressure  area  centered  over  the  Great  Lakes  caused  air  to  flow  down  the 
east  coast  to  the  Wallops  Island  site.  The  average  of  three  size  distributions  taken  between  1515  and 
1645  EDT  on  25  June  (centered  on  the  time  indicated  by  arrow  1  in  Fig.  16)  is  shown  as  curve  I  at  the 
left  in  Fig.  19,  The  total  number  of  particles  exceeded  15,000/cmV  During  the  next  three  days  the 
high-pressure  center  moved  eastward  to  a  position  off  the  New  England  coast,  as  shown  in  Fig.  17.  The 
air  arriving  at  the  measuring  site  had  increasingly  longer  trajectories  over  water  (and  most  likely 
different  source  regions).  The  change  in  the  size  distribution  during  this  period  is  shown  at  the  left  in 
Fig.  19,  at  times  which  correspond  to  the  meteorological  conditions  indicated  in  Fig,  16  by  the  arrows 
along  the  top  border.  By  1800  EDT  on  25  June,  the  local  wind  had  switched  from  north  to  northeas¬ 
terly,  and  the  total  count  had  dropped  from  15,200/cm’  to  12,700/cm’.  The  decrease  occurred  for  par¬ 
ticles  with  radii  smaller  than  0,03  ^m.  The  following  day,  26  June,  the  total  count  had  dropped  below 
1000/cm’  and  remained  between  800/cm’  and  1000/cm’  for  the  remainder  of  27  and  28  June.  The 
total  count  did  not  go  below  800/cm’.  The  size  distributions  shown  in  Fig.  19  are  averages  of  three 
individual  size  distributions  taken  during  a  90-minute  time  interval  centered  on  the  times  shown  in  Fig. 
16  except  for  the  size  distribution  for  1800  on  25  June  (arrow  2  in  Fig.  16),  which  was  an  average  of 
only  two  distributions.  The  dip  in  the  27  June  size  distribution  (curve  4)  at  about  0.06  fim  was  real 
and  existed  for  at  least  4  hours  but  was  not  a  persistent  feature  of  the  size  distribution  throughout  the 
low-count  interval  on  27  and  28  June.  When  this  dip  disappeared,  the  count  was  more  typically  about 
1000/cm’.  Also,  by  27  June  the  persistent  east  wind  resulted  in  higher  seas,  and  surf-generated  parti¬ 
cles  caused  a  significant  increa.se  in  the  number  of  particles  larger  than  about  1  pm. 
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Figure  19  shows  at  the  right  the  area  distributions  and  illustrates  the  advantages  gained  by  plotting 
diiTerent  moments  of  the  size  distribution.  The  area  distribution  emphasizes  the  differences  in  the  dis¬ 
tribution  at  large  sizes  and  has  a  tendency  to  decrease  the  steep  slopes  in  the  size  distribution  at  larger 
sizes. 

Unfortunately,  the  nephelometer  failed  on  25  June,  and  scattering  coefficients  were  not  recorded 
for  the  rest  of  the  deployment.  Some  spot  measurements  were  made  with  a  second  nephelometer  used 
in  the  variable-humidity  nephelometer  system. 

The  high  pressure  previously  described  brought  low  relative  humidities  and  good  visibilities,  and 
even  though  the  relative  humidity  increased  during  the  interval,  as  the  air  spent  a  longer  time  over 
water,  the  visibility  remained  high.  The  large  number  of  small  particles  (presumably  from  the  indus¬ 
trial  northeast)  recorded  on  25  June  had  little  effect  on  visibility.  Lower  visibility  (65  ~  2  x  lO'Vm) 
occurred  on  22  June.  The  weather  maps  indicated  that  this  air  had  passed  over  the  Carolinas  and 
southern  Virginia  before  arriving  at  the  site.  A  cluster  of  eight  size  and  area  distributions  taken  in  late 
afternoon  and  early  evening  on  22  June  are  shown  in  Fig.  20.  A  predominant,  relatively  stable  peak  in 
the  area  distribution  occurred  at  about  0.1  /am.  This  peak  radius  was  considerably  larger  than  the  peak 
radius  on  June  25  (Fig.  19).  The  area  distributions  have  not  only  a  predominant  maximum  but  a 
secondary  maximum  at  about  1  /am.  This  secondary  maximum  cannot  be  attributed  to  locally  generated 
surf  particles.  The  local  wind  at  the  time  was  from  the  southwest  and  could  not  have  transported  surf 
particles  from  the  beach,  which  was  80  m  east  of  the  site. 

Figure  21  shows  three  average  distributions  of  size  and  of  area  taken  on  21,  19,  and  22  June 
(curves  1,  2,  and  3).  The  scattering  coefficients  associated  with  these  three  times  were  0.40,  1.30,  and 
2.1  X  10  Vm  respectively.  The  values  of  the  extinction  coefficient  at  a  wavelength  of  0.55  f/m  calcu¬ 
lated  from  these  size  distributions  are  0.37.  1.42,  and  3.0  x  10  Vm  respectively.  As  would  be 
expected,  the  area  associated  with  the  number  of  particles  greater  than  0.1  /im  correlates  quite  well  with 
the  value  of  the  scattering  coefficient.  In  general,  we  observe  no  correlation  between  the  total  particle 
count  and  scattering,  because  the  total  count  (over  land)  is  dominated  by  the  smaller  particles.  The  air 
reaching  the  site  on  19  June  circulated  down  the  cast  coast  around  a  high-pressure  system  centered 
over  Ontario  and  probably  originated  over  Canada. 
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Shipboard  Measurements  off  the  East  Coast  of  the  United  States  (1980  USNS  Hayes  Cruise) 

In  April  1980  a  joint  cruise  of  the  NRL  Marine  Biology  Branch  and  the  NRL  Atmosplieric  Physics 
Branch  was  undertaken  off  the  east  coast  of  the  United  States.  The  Marine  Biology  Branch's  objectives 
required  stops  4  to  10  hours  in  duration  at  certain  stations  inside  and  outside  the  Gulf  Stream.  The 
purpose  of  the  Atmospheric  Physics  Branch’s  participation  was  to  obtain  data  on  aerosol  size  distribu¬ 
tions  along  the  east  coast  with  special  emphasis  on  the  evolution  of  the  size  distribution  in  air  advecting 
off  the  east  coast.  Our  objective  required  the  ship  to  be  heading  such  that  the  wind  speed  and  direction 
relative  to  the  ship  would  assure  uncontaminated  air  at  the  sampling  location. 
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The  ISNS  Haves  (I'-ACiOR  IM  is  a  ivsin-hullcil  (catamarani  ^ssssl  designed  es|ieeiall\  Im 
oeeaiK'graphie  roseareh.  The  portable  instrument  shelter  desenbed  eailiei  «as  installed  on  the  Hung 
bridge  of  the  ship  abo\e  (he  pilot  house  'l  ig  ddb  I  he  llni.g  bridge  lOMif  ol  the  pilot  luaisei  is  ahou: 
14  Ill  aboie  the  water.  The  air  sample  lor  the  mobiliti  anab/ei  aiui  (  N  counter  was  taken  ihinagh  a 
stainless  steel  tube  shaped  like  a  shepherrl's  stall' and  mounted  with  the  mlet  about  '  ni  abioe  the  Iront 
of  the  lulot  house.  The  Roieo  optieal  sensor  was  on  the  leailing  edge  o!  the  metal  stand  on  the  liorii  o| 
the  shelter  .Airilow'  measurements  on  a  model  of  the  I'SNS  //uus  were  (irenousli  made  m  a  wind 
tunnel  at  the  Nasal  Ship  Research  and  Development  (enter,  Hethesda.  Maiiland  iJeik.  14’' i  figure 
23  shows  the  results  of  these  tests  and  indicates  that  the  sampling  sites  3  ’o  3  m  abose  the  leading 
edge  of  the  pilot  house,  provided  good  e.xposure  for  sampling  when  the  lelaiisc  wind  was  trom  the  for¬ 
ward  sector  In  practice  the  CN  counter  provides  a  sensitise  test  of  ship's  contamination  in  the  clean 
ocean  environment  F.ven  small  traces  of  ship's  pollution  (from  exhaust  or  g.illcs  sentst  produce  large 
numbers  of  CN,  easily  detected  by  the  CN  counter  No  indication  of  ship's  pollution  was  eser  detected 
unless  the  relative  wind  was  more  than  70°  to  either  side  of  the  ship's  heading 

The  cruise  began  on  the  afternoon  of  10  -Xpril  148t)  and  followed  the  track  shown  n  fig  24  I  he 
numbered  boxes  are  stations  where  the  ship  stopped  for  2  to  12  hours  for  the  obseiw .iiions  b\  the 
marine  biologists.  Station  2  was  a  reference  station  visited  on  three  occ.ismns.  and  station  3  was  visited 
twice.  Contamination  originating  from  the  ship  itself  usuallv  made  it  impossible  to  sample  when  the 
ship  was  on  station. 

Time  plots  of  radon,  CN,  CCN,  the  scattering  coelTtcient.  .md  meteorological  d,iia  .ue  shown  m 
Fig.  25a  through  25d.  Data  for  15  .April  are  omitted  because  most  of  the  Cciuipmcni  w.ic  not  oper,iting 
on  that  day.  Most  of  the  breaks  in  the  data  occurred  when  the  shi  '  w,is  on  station 
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After  the  ship  left  port  on  10  April  and  continuing  through  12  April  the  cruise  area  was  under  the 
influence  of  a  well-defined  but  rather  weak  (5-knot)  flow  from  off  the  Virginia  and/or  North  Carolina 
coast.  Air-mass  trajectories  for  this  interval  are  shown  in  Fig.  26.  The  relatively  uniform  radon  con¬ 
centrations  throughout  this  interval  (Fig.  25a)  attest  to  the  uniformity  of  the  air  mass  on  a  synoptic 
scale. 


l  ig.  2(i  —  I  rjiccliiric^  showing  Ihc  origin  oC  ihc  .iir  in  which  the  tjisuihii- 
lions  wtYt'  mcjsurcil  >il  ID)  on  )l).  Jl,  12.  iinJ  1.2  .ypril  .il  the  loc.i- 
lions  indic.ilctl  h\  (illciJ  circles 


A  well-defined  cold  front  moved  through  the  operations  area  during  the  night  of  13-14  April.  The 
radon,  scattering-coefficient,  and  CCN  data  patterns  are  similar  to  those  during  previous  transits  of 
frontal  areas  over  the  ocean  (Larson  and  Bressan,  1980).  The  front  is  shown  as  a  line  on  weather 
maps,  indicating  the  location  of  the  change  in  wind  direction,  but  radon  data  suggest  that  the  frontal 
area  is  a  broad  zone  where  the  continental  materials  are  greatly  enhanced  as  far  as  a  few  hundred 
kilometers  ahead  of  the  front.  Radon,  aerosol  scattering,  CCN,  and  other  indicators  of  continental  air 
exhibit  much  fine  structure  in  these  frontal  zones.  Sometimes  the  fine  structure  is  very  similar  from 
one  variable  to  the  next,  but  in  some  cases,  such  as  this,  there  is  little  correlation  other  than  the  overall 
continental  enhancement.  The  increase  in  radon  concentration  appreciably  lags  the  sharp  increases  in 
aerosol  scattering  and  CCN  and  CN  concentrations  starting  about  1500  EDT  on  13  April,  and  except 
for  a  small  change  in  scattering  the  abrupt  change  in  wind  direction  marking  the  frontal  boundary  is  not 
accompanied  by  changes  in  other  indicators  of  continental  air. 

By  0900  EDT,  14  April,  the  front  has  changed  to  a  warm  front,  and  the  change  in  wind  direction 
from  east  to  south  is  a  consequence  of  the  approach  of  another  low-pressure  area  over  West  Virginia. 
Between  1500  and  1800  EDT  on  14  April,  the  scattering  coefficient  and  concentrations  of  radon  and 
particulates  dropped  to  low  values,  and  the  ship  was  in  clean  maritime  air  brought  into  the  area  by  the 
combined  influence  of  the  low  over  the  eastern  United  States  and  a  high  pressure  east  of  Bermuda. 
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Data  for  15  April  are  not  included  because  most  of  the  instrumentation  was  shut  down  due  to 
heavy  seas.  Radon  data  showed  maritime  air  until  about  0900  EDT.  During  the  remainder  of  the  day 
the  radon  concentrations  were  in  the  range  It)  to  15  pCi/m’  that  prevailed  through  noon  on  16  April. 
A  cold  front  from  the  west  passed  the  ship  late  on  15  April,  and  another  cold  front  accompanied  by 
light  rain  passed  over  the  ship  just  before  arrival  at  station  II  at  2000  EDT  on  16  April.  The  peak  in 
total  aerosol  count  and  in  CCN  on  16  April,  while  scattering  remained  relatively  uniform,  can  be  attri¬ 
buted  to  larger  numbers  of  particles  with  radii  less  than  0.07  jim  which  do  not  contribute  to  the  scatter¬ 
ing. 


After  passage  of  the  front  on  16  April,  the  ship  came  under  the  influence  of  a  high-pressure  sys¬ 
tem  The  center  of  the  high  pressure  moved  from  over  the  Great  Lakes  to  just  off  Cape  Hatteras  on  17 
.April  and  remained  stationary  for  the  next  3  days.  The  combination  of  brisk  northerly  winds  and  rela¬ 
tively  low  radon  and  total  particles  is  consistent  with  an  assessment  based  on  the  weather  maps  that  the 
air  reaching  the  ship  on  the  morning  of  the  17  April  was  originally  maritime  air  brought  into  the  area 
after  a  quick  transit  over  New  England  by  the  high  over  the  Cireat  Lakes  and  a  low  off  Newfoundland 
The  fast  movement  minimizes  the  accumulation  of  continental  additives  to  the  air.  The  winds  became 
more  westerly  in  the  afternoon,  suggesting  that  the  large  increase  in  small  particles  (increased  total 
count  but  relatively  unchanged  scattering)  results  from  industrial  pollution  from  the  United  States  or 
Canada.  The  radon  concentrations  were  relatively  low  and  variable  on  18,  19,  and  20  April,  when  the 
ship  was  in  subsiding  air  near  the  center  of  the  high  pressure. 

•A  cold  front  crossed  the  ship's  location  about  noon  on  21  April,  accompanied  b)’  a  change  in  wind 
direction  and  a  sharp  increase  in  radon,  indicating  that  the  air  behind  the  Iront  was  continental.  The 
ship  remained  in  mostly  continental  air  throughout  the  remainder  of  the  cruise,  near  the  coast  and  in 
the  Chesapeake  Bay. 

Air  trajectories  were  computed  using  the  NOAA  regional-continental  scale  transport,  diffusion, 
and  deposition  model  (Heffter  et  al,,  1975).  This  model  has  worked  well  over  land,  where  sufficient 
data  are  available,  and  it  has  worked  reasonably  well  in  remote  ocean  areas  in  that  it  confirmed  conclu¬ 
sions  about  the  continental  or  maritime  character  of  air  based  on  measured  radon  concentrations 
Because  the  model  does  not  take  into  account  vertical  air  motions  or  discontinuities  at  fronts,  and 
because  data  are  few  for  areas  over  the  ocean,  the  near-coastal  area  of  this  cruise,  with  the  resulting 
importance  of  mesoscale  air  movements,  did  not  suit  itself  well  for  use  of  this  model.  However,  in 
several  cases  the  computed  trajectories  confirmed  the  conclusions  that  could  be  drawn  from  particle 
data,  radon  data,  and  weather  maps  as  to  the  origin  and  history  of  the  air  at  the  sample  collection  site. 
The  trajectories  from  this  model  are  used  in  this  report  only  when  they  are  supported  by  other 
meteorological  data. 

Overview  oj  the  Size  Distribution  Measurements 

Over  300  aerosol  size  distributions  were  obtained  during  the  cruise.  An  analysis  of  each  size  dis- 
trihution  in  conjunction  with  the  meteorological  history  of  the  air  mass  is  not  possible.  In  this  subsec¬ 
tion  some  gener.il  characteristics  of  the  size  distributions  obtained  during  the  cruise  will  bo  described 
and  in  the  following  subsections  case  studies  of  several  interesting  intervals  will  be  prescnied. 

The  size  distribution  over  the  ocean  was  often  observed  to  remtiin  remarkablv  stable  with  little 
change  observed  over  the  interval  of  hours.  This  is  in  contrast  to  size  distributions  over  land,  where 
temporai  variations  over  similar  time  intervals  arc  quite  large,  f'igure  27  is  a  composite  of  five  consecu¬ 
tive  size  and  area  distributions  taken  during  2-1/2  hours  (0430  to  0700  EDT  on  11  April)  just  prior  to 
arrival  on  station  2.  The  five  distributions  are  for  the  most  part  indistinguishable.  .At  other  times  the 
distributions  would  change  considerably  over  an  interval  of  hours.  These  changes  were  usually ,  but  not 
iilways,  associated  with  some  observable  meteorological  changes,  with  the  distance  the  air  traveled  over 
water  before  arriving  at  the  ship,  or  with  the  .sea  state  Since  the  distributions  often  remained  stable. 
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1  ij:  27  —  ( dnipiiMio  ot  ti\c  cunsccuiivc  m/c  vlisu i^uiinns  .tnd  ‘d  ihc  concspcndiny  .irca  disirihuiu'ns 
taken  during  a  2-l/2-h»UJr  interval.  illuNtraiine  the  stability  uf  the  m/c  djsjnhuii«>n  oflen  observed  over 
the  ocean 


the  data  were  sorted  into  periods  when  little  or  no  change  occurred,  and  average  distributions  for  these 
periods  were  calculated.  This  was  an  effective  way  of  reducing  the  numbers  of  data  and  of  increasing 
their  statistical  significance.  In  the  next  subsections  the  distributions  presented  will  be  averages  of  two 
or  more  distributions. 

The  variability  in  the  size  and  area  distributions  over  longer  intervals  is  illustrated  in  Figs.  28.  29. 
and  30  where  each  curve  is  an  individual  size  or  area  distribution.  The  left-hand  plot  in  Fig.  28  is  a 
composite  of  75  size  distributions  taken  between  the  time  the  ship  left  the  Chesapeake  Bay  (1700  EDT 
on  10  April)  and  0100  EDT  on  14  April,  at  a  point  between  stations  6  and  7.  The  variability  in  the 
small  particles  covers  about  four  orders  of  magnitude.  The  variability  in  the  size  distribution  is  remark¬ 
ably  low  between  0.07  and  0.5  fim  in  radius.  The  region  between  about  0.25  a. id  0.4  ^im  is  the  region 
of  overlap  between  the  mobility  analyzer  and  the  optical  particle  counter,  and  the  spread  in  this  region 
is  artificially  increased  because  of  the  difference  in  the  size  distribution  measured  by  the  two  instru¬ 
ments.  When  the  value  from  the  largest-particle  channel  is  below  the  abscissa,  it  is  not  plotted  This  is 
the  reason  the  size  distribution  sometimes  stops  at  1.2  jam  rather  than  about  2  /um.  The  average  slope 
for  the  size  range  between  0.1  and  1  /im  is  nearly  -4  during  this  time,  in  agreement  with  the  slope  ori¬ 
ginally  proposed  by  Junge  (1963). 

The  steep  slope  in  the  size  distribution  at  radii  greater  than  0.08  /am  tends  to  obscure  the  varia¬ 
tions  in  numbers  at  the  larger  radii.  At  the  right  in  Fig.  28  a  composite  of  the  same  data  plotted  as  area 
distributions  is  presented.  This  plot  decreases  the  slope  and  magnifies  the  variations  at  larger  sizes.  It 
also  magnifies  the  mismatch  in  the  distributions  obtained  with  the  mobility  analyzer  and  the  optical  par¬ 
ticle  counter  in  the  overlap  region  between  about  0.2  /um  and  0.4  /um. 

The  air  during  this  period  was  from  the  continental  United  States  but  with  differing  trajectories 
over  water.  Figure  26  shows  typical  trajectories  during  this  portion  of  the  cruise.  The  size  distributions 
with  the  largest  numbers  of  small  particles  were  taken  shortly  after  the  ship  had  left  Chesapeake  Bay  in 
air  that  had  a  short  trajectory  over  water.  As  will  be  discu.ssed  later,  there  is  apparently  always  an  initial 
rapid  decrease  in  the  number  of  small  particles  as  the  air  advccts  from  land  out  over  the  ocean.  How¬ 
ever,  beyond  a  certain  point  it  is  not  clear  that  the  decrease  is  a  simple  function  of  distance  (or  time) 
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l  ii;  30  —  Composite  of  104  dislnhutions  taken  between  1'^  April  and  31  Apiil 

over  water.  (Part  of  the  problem  here  is  the  difficulty  of  establishing  reliable  back-trajectories  of  sur¬ 
face  air  at  the  ship  for  longer  times  and  the  difference  in  sources  of  the  different  samples  in  that  the 
measurements  are  not  on  the  same  trajectory.) 

The  air  encountered  during  this  portion  of  the  cruise  had  apparently  not  gone  through  a  precipita¬ 
tion  cycle  during  its  transit  from  land  to  the  ship.  This  probably  plays  an  important  role  in  explaining 
the  invariance  of  the  size  distribution  between  0.07  fxm  and  0.5  and  the  large  number  of  particles 
in  the  size  range  O.OS  to  O.I  /utn 

Figure  29  shows  a  composite  of  40  size  distributions  taken  between  0200  EDT  on  14  April  and 
1400  EDT  on  16  April  (no  data  were  taken  on  15  April  and  on  the  morning  of  16  April  because  of  a 
storm).  Again  there  is  a  large  variation  in  small  particles.  Here  the  stability  observed  in  the  range  0.07 
to  0.5  /am  in  Fig.  28  has  degenerated.  The  size  distributions  v^ith  low  concentrations  in  this  intermedi¬ 
ate  size  range  were  probably  the  result  of  wet  removal  processes  associated  with  a  frontal  system 
These  same  size  distributions  exhibit  a  large  number  of  particles  greater  than  1  ^/m  These  larger  parti¬ 
cles  are  surface-generated  by  high  seas. 

Figure  30  is  a  composite  of  104  size  distributions  taken  during  the  latter  part  of  the  cruise 
between  0800  on  17  April  and  1300  on  21  April.  This  interval  was  dominated  by  a  nearly  stationary 
high-pressure  system  moving  slowly  off  the  east  coast.  Again  there  is  a  large  variability  in  the  small 
particle  concentrations,  and  the  slope  between  0.1  /xm  and  1  /am  is  about  -4. 

Evolution  of  the  Size  Distribution  in  Air  Masses  Advecting  from  Land  to  Sea 

One  important  objective  of  our  participation  in  the  1980  cruise  was  to  observe  the  change  in  the 
size  distribution  of  small  particles  in  air  masses  as  they  advect  from  the  east  coast  out  over  the  Atlantic. 
The  initial  decrease  in  the  small  aerosol  population  as  continental  air  advects  out  over  the  water  is 
clearly  evident  in  a  number  of  cases  from  our  data.  However,  determination  of  the  precise  rate  of 
decay  along  a  given  trajectory  as  a  function  of  particle  size  is  not  possible  from  the  data  obtained.  The 
main  difficulties  are  that  the  complex  meteorological  processes  cause  uncertainties  in  the  accuracy  of 
the  calculated  back-trajectories  of  air  encountered  at  the  surface  and  that  measurements  were  not  along 
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a  single  trajectory.  In  many  cases  the  .calculated  back-trajectories  (for  an  altitude  of  450  m)  did  not 
agree  with  analysis  of  surface  weather  charts  or  local  observation  of  wind  and  r^don  concentrations 
Despite  the  lack  of  knowledge  of  the  exact  trajectories  of  the  surface  air,  interesting  data  on  the  deca> 
were  obtained.  Two  case  studies  will  be  presented:  between  port  and  station  2  and  between  stations  ' 
and  9. 


The  ship  left  Norfolk  harbor  about  1400  EDT  on  10  April,  passed  through  the  entrance  to  the 
Chesapeake  Bay  at  about  1530  and  headed  east  for  station  2  (Fig.  24).  The  night  before,  a  front 
accompanied  by  rain  passed  through  the  area  and  out  to  sea.  The  air  on  the  morning  of  10  April  was 
dry  and  clear.  After  the  ship  left  the  Bay,  the  air  temperature  dropped  during  the  first  5  hours  from 
27°C  to  13°C  and  then  increased  to  about  21°C  as  the  ship  approached  station  2.  The  latter  warming 
was  attributed  to  the  influence  of  the  Gulf  Stream.  Curve  1  in  the  left-hand  plot  of  Fig.  31  shows  the 
size  distribution  taken  in  the  Chesapeake  Bay  between  the  James  River  and  the  mouth  of  the  Bay.  The 
total  particulate  count  was  about  30,000/cm^  and  highly  variable. 
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Curve  2  is  the  average  of  ten  size  distributions  taken  between  1700  FDT  and  2100  FDT  about  15 
to  60  nmi  off  shore.  The  total  particle  count  was  about  9000/cmT  During  this  interval  there  was  little 
systematic  change  in  the  size  distribution.  Unfortunately,  the  ship’s  speed  was  not  recorded  during  this 
interval,  so  true  wind  speed  and  direction  could  not  be  accurately  calculated.  From  estimates  of  the 
ship’s  speed  the  estimated  wind  direction  was  from  the  south-southwest.  As  the  ship  moved  out  to  sea, 
the  air  encountered  undoubtedly  had  a  back-trajectory  for  which  landfall  was  progressively  to  the  south. 
For  curves  3  and  4  landfall  was  probably  somewhere  along  the  North  Carolina  coast. 

The  average  of  three  distributions  taken  between  2330  and  0100  FDT  on  1 1  April  about  60  to  90 
nmi  off  the  coast  is  shown  as  curve  3.  Curve  4  is  the  average  of  five  distributions  taken  between  0430 
and  0700  FDT  near  station  2,  about  130  nmi  ofT  the  coast. 

Figure  32  shows  the  calculated  air-mass  trajectories  supplied  by  NOAA  for  450-m  level  during  this 
interval,  (Surface  winds  would  indicate  that  the  surface  air  came  more  from  the  south  than  indicated 
by  the  trajectories  for  450  m.)  The  weather  was  clear,  and  there  was  no  indication  that  the  air  had 
passed  through  any  region  of  precipitation  during  the  previous  24  hours,  and  certainly  not  through  any 
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precipitation  after  the  air  had  been  over  water.  If  these  trajectories  are  used  to  calculate  the  time  the 
air  spent  over  water  between  the  coast  and  arrival  at  the  ship,  the  values  are  15  to  3  hours  for  curve  2. 
5  hours  for  curve  3,  and  9.5  hours  for  curve  4.  It  is  suspected  that  these  limes  calculated  I'rnm  the 
450-m-level  trajectories  are  minimal  and  that  the  real  times  are  somewhat  longer 

The  decay  in  the  number  of  small  particles  as  the  air  spends  more  time  over  water  is  obvious 
from  Fig.  31  and  is  undoubtedly  due  to  the  loss  of  small  particles  by  coagulation  with  each  other  and 
with  larger  particles  and  by  growth  of  smaller  particles  to  larger  sizes  by  condensation  of  gas-phase  reac¬ 
tion  products  on  their  surfaces  (details  of  the  latter  mechanism  being  given  in  Hoppel.  1975)  These 
mechanisms  have  little  effect  on  the  particles  greater  than  0.08  yam  in  radius.  The  number  of  large  par¬ 
ticles  first  increased  as  the  ship  moved  into  the  open  sea  and  then  decreased  continually  until  the  ship 
reached  station  2.  This  decrease  in  large  particles  between  curves  2  and  4  is  more  difficult  to  under¬ 
stand.  The  decrease  occurred  while  the  wind  also  decreased  from  about  16  knots  to  7  knots,  which  sug¬ 
gests  the  large  particles  may  be  wind-generated  particles  originating  at  the  sea  surface.  However,  this 
explanation  is  not  entirely  satisfactory,  because  at  other  times  on  the  cruise  wind-speed  changes  of  this 
magnitude  did  not  cause  a  large  change  in  the  larger  particles.  The  change  is  probablv  related  not  only 
to  the  instantaneous  wind  speed  but  also  to  the  wind  and  sea-state  history  over  the  prior  24  hours  or  so 
(Gathman,  1983). 

The  measured  scattering  coefficients  at  the  times  curves  2,  3,  and  4  were  obtained  were  0  56, 
0.47,  and  0.58  x  10  Vm,  respectively.  This  would  indicate  that  the  scattering  (at  0  55  /urn  in 
wavelength)  was  not  greatly  affected  by  the  particles  larger  than  0.5  )am  in  radius.  This  is  in  agreement 
with  the  detailed  analysis  of  extinction  discussed  elsewhere  in  this  report  Particles  smaller  than  0  3  /jm 
are  predominantly  responsible  for  the  scattering  of  visible  light. 

Another  particularly  interesting  period  of  the  cruise  occurred  between  stations  5  and  9  (Fig.  24) 
from  1500  EDT  on  13  April  to  1100  EDT  on  15  April.  The  number  of  small  particles  increased  as  the 
ship  approached  the  coast  and  then  decreased  again  as  the  ship  sailed  away  from  the  coast  As  the  ship 
left  station  8,  a  storm  was  encountered,  and  the  aerosol  size  distribution  associated  with  the  storm  was 
measured 
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During  this  2-da\  interval  the  meteorology  was  fairly  complex.  At  the  start  of  the  interval  a  cold 
front,  trailing  from  a  low  over  Nova  Scotia,  passed  across  the  ship's  position.  This  front  connected  to 
another  low-pressure  system  centered  over  southern  Mississippi.  A  ridge  of  high  pressure  separated  the 
two  lows.  During  the  transit  from  station  5  to  station  7  the  ship  was  in  air  behind  the  cold  front.  This 
air  was  presumably  of  recent  continental  origin.  The  true  winds  were  too  light  to  calculate  accuratelv 
during  most  of  the  transit  to  shore  but  were  generally  from  the  west.  The  low  that  was  over  the  .Missis¬ 
sippi  Gulf  Coast  deepened  somewhat  as  it  tracked  to  the  north-northeast  and  was  centered  over  western 
Kentucky  by  0800  EDT  on  14  April.  The  circulation  around  this  low  caused  the  wind  at  the  ship  to 
strengthen  and  shift  to  the  east  during  the  early  morning  hours  of  14  April.  When  the  ship  arrived  at 
station  7,  the  winds  were  easterly  at  22  knots.  The  skies  had  become  increasingly  overcast  during  the 
shoreward  part  of  the  track,  with  heavy  overcast  and  some  mist  where  the  track  came  closest  to  the 
coast. 


Figure  33  shows  four  average  size  distributions  and  area  distributions  taken  between  stations  5 
and  7  and  clearly  shows  the  increase  in  small  particles  as  the  ship  got  closer  to  land.  Curve  1  in  each 
plot  is  the  average  of  four  distributions  taken  between  2000  and  2230  on  13  April  when  the  ship  was 
between  stations  5  and  6,  The  total  count  was  about  3000/cm^.  Curve  2  in  each  plot  is  the  average  of 
six  size  distributions  taken  between  0350  and  0500  on  14  April.  The  total  count  was  about  3500/cm' 
Curve  3  in  each  plot  is  the  average  of  two  distributions  between  0600  and  0700  at  the  closest  point  to 
land.  The  total  count  was  again  about  3500/cm',  with  the  increase  in  small  particles  being  balanced  by 
the  decrease  in  the  size  range  0.05  to  0.2  /xm.  Curve  4  in  each  plot  is  the  average  of  two  distributions 
taken  between  0700  and  0845  just  prior  to  the  ship’s  arrival  at  station  7.  The  total  count  was  about 
5000/cm'.  When  the  data  for  curves  3  and  4  were  obtained,  the  wind  had  intensified  to  about  22  knots 
and  was  shifting  from  northwest  to  easterly. 
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The  ship  remained  at  station  7  from  about  0800  to  1400  EDT,  during  which  time  no  data  W'ere 
taken.  At  1400  the  ship  left  station  7  and  slow"!)’  transited  to  station  8,  where  it  remained  from  1800  to 
2400  EDT.  Throughout  this  interval  the  seas  became  rougher,  with  the  wind  speed  ranging  from  17  to 

25  knots  from  the  south-southeast.  As  the  ship  sailed  from  station  7  to  station  8,  the  air  became 

noticeably  cleaner,  as  indicated  by  the  lower  particle  counts,  lower  radon  concentrations,  and  smaller 
values  of  the  scattering  coefficient.  The  surface  weather  map  showed  that  the  warm  front  associated 
with  the  low  over  the  eastern  United  States  moved  north  across  the  ship's  position  at  about  1430  EDT 
on  14  April.  Curves  1  of  Fig.  34  are  the  same  as  curves  4  of  Fig.  33  and  are  the  distributions  upon  the 
ship's  arrival  at  station  7.  Curve  2  in  each  plot  is  the  average  of  two  distributions  taken  between  1630 
and  1730,  just  prior  to  the  ship's  arrival  at  station  8.  The  total  particulate  count  had  dropped  to 

1400/cm^,  but  there  was  an  increase  in  particles  larger  than  1  /urn  in  radius,  believed  to  be  particles 

generated  by  the  higher  seas. 

When  the  ship  left  station  8  at  0000  EDT  on  15  April,  the  winds  were  very  strong  (30  to  35 
knots)  from  the  southeast.  The  seas  became  increasingly  rough,  and  the  wind  became  so  strong  that 
there  was  a  reverse  flow  of  air  through  the  nephelometer  due  to  Bernoulli  suction  of  the  wind  over  the 
air  inlet.  The  relative  humidity  had  increased  and  was  between  95'’/(i  and  lOO'/n.  Curve  3  in  each  plot  is 
the  average  of  two  distributions  taken  during  this  interval  (0240  to  0330  EDT).  The  total  particle  count 
dropped  to  350/cm^,  the  lowest  recorded  on  the  cruise.  There  was  a  further  increase  in  the  number  of 
particles  larger  than  0.5  ^im  in  radius.  These  larger  particles  are  believed  to  be  salt  particles  generated 
by  the  high  seas  and  wind.  Rain  and  fog  developed  shortly  after  the  curve-3  data  were  obtained. 
Because  of  increasingly  stormy  conditions,  the  instruments  were  turned  off  at  about  0400,  anu  all 
efforts  were  spent  securing  and  protecting  the  equipment  from  damage.  The  presence  of  the  shelter  on 
the  flying  bridge  makes  for  a  particularly  rough  ride  in  heavy  seas  because  of  the  long  lever  arm 
between  the  ship's  center  of  gravity  and  the  bridge.  Although  the  storm  passed  that  night,  the  seas 
remained  heavy,  and  no  data  were  taken  for  the  next  1-1/2  days. 

The  dramatic  change  in  the  size  distribution  between  curves  1  and  2  in  Fig,  34  is  associated  with 
the  change  in  air  mass  following  the  passage  of  the  warm  front.  We  believed  that  prior  to  the  frontal 
passage  the  instruments  were  sampling  air  which  recently  had  been  over  the  continent.  With  the  pas¬ 
sage  of  the  warm  front  this  air  was  replaced  by  air  which  had  a  longer  trajectory  over  water,  as  can  be 
seen  from  Fig.  35,  which  shows  the  NO  A  A  air-mass  trajectories  for  2000  EDT  on  14  April  and  for  0200 
on  15  April.  These  trajectories  are  consistent  with  the  synoptic  situation  which  indicated  that  following 
the  passage  of  the  warm  front  the  air  reaching  the  ship  resulted  from  flow  on  the  back  side  of  a  broad 
high-pressure  area  (102,8  kPa  or  1028  mb)  with  one  center  near  30°N,  55°W.  The  low  count  in  the 
third  size  distribution  probably  indicates  air  which  underwent  a  wet-scavenging  cycle.  Radon  concentra¬ 
tions  of  2  pCi/m^  indicate  that  the  air  had  not  been  over  land  for  some  time. 

Figures  36  and  37  are  plots  of  the  extinction  coefficients  computed  for  the  size  distributions  given 
by  curves  1  and  3  in  the  left-hand  plot  of  Fig.  34.  The  dramatic  change  in  the  size  distribution  is 
reflected  in  the  relative  contribution  of  the  particles  in  each  size  range  to  the  total  extinction.  Curve  3 
in  the  left-hand  plot  of  Fig.  34  is  the  only  size  distribution  we  encountered  such  that  particles  in  the 
size  range  0.9  to  3  ^m  made  a  major  contribution  (in  this  case  the  largest  contribution)  to  extinction  at 
0.55  /am  in  wavelength  while  particles  smaller  than  about  0.3  /am  made  a  negligible  contribution. 

Size  Distributions  Under  High  Pressure  Off  the  East  Coast  of  the  United  States 

After  the  storm  the  seas  remained  rough  on  15  April,  with  winds  from  the  south  at  about  26 
knots.  Late  in  the  day  on  15  April,  a  cold  front  from  the  west  passed  over  the  ship,  and  the  winds 
shifted  more  to  the  west  at  about  22  to  30  knots.  The  equipment  which  had  been  turned  off  early  on 
15  April  during  the  storm  was  turned  on  again  while  the  ship  was  on  station  4.  Another  cold  front 
accompanied  by  light  rain  passed  over  the  ship  just  before  its  arrival  on  station  1 1  at  2000  EDT  on  16 
April  For  the  next  4-1/2  days  the  ship  was  under  a  high-pressure  system  which  had  moved  down  from 
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the  southern  Great  Lakes  region  on  16  April  to  over  Cape  Halteras  on  17  April  At  the  coast,  the 
movement  slowed,  and  the  high  moved  slowlv  out  to  sea,  degenerating  into  a  high-pressure  ridge  b>  19 
April. 


Curve  1  in  each  plot  of  Fig,  38  is  the  average  of  eight  distributions  taken  during  a  5-hour  interval 
before  the  ship  arrived  at  station  12.  Each  of  the  eight  distributions  had  the  peculiar  double  peak  exhi¬ 
bited  in  curve  1.  The  wind  during  this  interval  was  from  the  northwest  and  was  decreasing  in  speed 
from  about  26  knots  to  16  knots.  The  total  particulate  count  was  1500  to  2000/cm\  While  the  ship 
was  on  station  12,  the  total  particulate  count  increased,  and  it  remained  high  as  the  ship  approached 
closer  to  land.  Curve  2  in  each  plot  of  Fig.  38  is  the  average  of  eight  size  distributions  taken  between 
station  12  and  the  closest  approach  to  land.  The  total  particulate  count  here  was  about  7000  to 
10,000/cm'.  The  winds  were  from  the  west  at  about  9  knots,  and  the  ship  was  just  north  of  the  center 
of  the  high  pressure.  From  the  available  data  it  is  difficult  to  determine  the  origin  of  the  air.  The  tra¬ 
jectory  for  the  450-m  level  shows  air  flowing  south  from  the  North  Atlantic  and  coming  close  to,  but 
never  quite  over,  land.  The  trajectory  for  900  m  shows  the  air  coming  from  the  opposite  direction 
(from  the  Bahama  Islands).  The  local  winds  were  front  the  west-southwest.  Although  the  increase  in 
the  total  aerosol  count  would  indicate  a  recent  trajectory  over  land,  the  radon  concentration  was  quite 
low,  indicating  that  the  air  had  not  spent  much  time  over  land.  Our  best  estimate  is  that  the  surface  air 
had  come  in  from  the  North  Atlantic  and  had  a  short  trajectory  over  land  (or  had  mixed  with  continen¬ 
tal  air)  before  going  back  over  the  sea  to  the  location  of  the  ship.  The  air  may  further  have  been 
diluted  with  subsiding  air  under  the  high-pressure  system. 
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f  ig  38  —  Three  average  disiribuiions  (1)  average  of  eight  distributions  taken  lust  prior  to  the  ship's  arrival 
on  station  12,  with  all  eight  cshibiimg  the  double  peak,  (2)  average  of  eight  distributions  between  station  12 
and  ihc  eloscsi  approach  to  land,  and  (3)  distributions  on  the  ship's  arrival  at  station  2  at  0700  HDT  on  18 
April 


As  the  ship  turned  southeast  away  from  land  at  about  2300  EOT  on  17  April  and  headed  toward 
stations  2  and  13,  the  number  of  small  particles  began  a  slow  decrease.  When  the  ship  arrived  at  sta¬ 
tion  2,  7  hours  later,  the  si/>’  distribution  was  as  shown  by  curve  3  in  the  left-hand  plot  with  a  total  par¬ 
ticulate  count  of  about  .3000/cm',  The  air  was  dry  compared  to  ty  pical  oceanic  humidities,  reflecting 
the  influence  of  the  high  pressure.  The  rei  oivc  humidity  was  about  60'Vli  at  station  12  and  dropped  to 
37"/ii  on  station  2.  when  the  ship  was  at  the  center  of  the  high. 
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The  distributions  obtained  from  18  April  to  20  April  are  diffieult  to  interpret  in  terni->  of  the 
meteorology  and  air-mass  origin.  In  Fig.  30  the  average  distributions  during  two  periods  when  the  dis¬ 
tributions  remained  quite  stable  are  shown.  Curve  1  in  each  plot  is  the  average  of  seven  si/e  distribu¬ 
tions  taken  between  1700  and  2400  on  18  .April  near  and  on  station  13.  .All  of  the  si/e  distributions 
individually  exhibited  the  double  peak.  The  wind  was  weak  from  the  north-northwest,  seas  were  calm, 
and  the  particulate  count  was  4500  to  5500/emV 
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Curve  2  in  each  plot  of  Fig.  39  is  the  average  of  16  similar  distributions  measured  between  1430 
EOT  on  20  April  and  0030  on  21  April  (the  first  half  of  the  transit  between  stations  I.''  and  16)  The 
wind  was  from  the  south  at  about  13  knots,  and  the  total  count  was  in  the  range  1500  to  2000/cnT 
Mere  the  peak  of  the  area  distribution  occurred  at  a  radius  exceeding  0.1  /am.  which  was  about  the  larg¬ 
est  radius  at  which  the  peak  occurred  on  the  cruise  and  much  larger  than  ever  observed  over  land 

Return  to  the  Chesapeake  Bay 

The  weather  maps  for  21  April  showed  a  cold  front  moving  across  the  location  of  the  ship  in  early 
afternoon  while  the  ship  was  on  station  16,  about  40  nmi  off  Cape  Ilatteras  There  was  no  significant 
change  in  temperature  or  cloudiness  associated  with  the  passage  of  the  fiont  that  was  observed  at  the 
ship.  However,  there  was  a  significant  increase  in  wind  speed,  scattering  coefficient,  and  total  aerosol 
count  (which  increased  from  1500  to  5000/cm’)  as  shown  in  Fig.  25d  at  about  1300  EDT  on  April  21 
The  increase  in  radon  concentrations  from  about  3  to  25  pCi/m’  confirms  that  after  the  passage  of  the 
front  the  air  had  a  more  recent  (and/or  longer)  trajectory  over  land  than  air  preceding  the  front  (This 
change  was  not  evident  from  NOAA  air-mass  trajectories.)  The  change  in  the  size  and  area  distribu¬ 
tions  associated  with  the  passage  of  the  front  is  shown  in  Fig.  40,  where  curve  1  in  each  plot  is  the 
average  of  three  distrl'  .Aons  taker  between  0715  and  0830  EOT  just  prior  to  the  ship's  arrival  on  sta¬ 
tion  16  and  curve  2  is  me  average  of  six  distributions  taken  between  1330  and  1650  EOT  after  the  ship 
left  station  16  en  route  to  the  mouth  of  the  Bay. 

.Starting  at  1600  EDT  on  21  April  the  ship  sailed  eastward  away  from  the  C  hesapeake  Bay  until 
2300  EDT,  when  it  turned  180°  and  headed  back  toward  the  Bay.  During  most  of  this  interval,  the 
total  aerosol  count  remained  in  the  range  3000  to  5000/cm',  and  the  size  distribution  was  similar  to 
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that  shown  b\  curve  2  in  the  lel'l-hand  plot  of  l  ig,  40.  During  the  6  hours  from  0600  to  1200  iiDT  (22 
April)  as  the  ship  approached  the  mouth  of  the  Ba>  and  proceeded  to  station  17  inside  the  Ba>.  there 
was  a  dramatic  decrease  in  the  aerosol  scattering  coefficient  and  a  corresponding  change  in  the  si/e  dis¬ 
tribution.  The  number  of  particles  larger  than  0.06  /am  decreased  dramaticalfs .  and  the  number  of 
small  particles  increased,  such  that  the  total  concentration  increased  from  about  6000/cm ‘  to  about 
20,000/cm'.  The  temperature  increased  during  this  time,  and  the  relative  humiditv  decreased  from 
about  70'/,  to  about  50'v:,.  There  were  no  navigational  data  available  from  which  to  calculate  wind  speed 
and  direction  during  transit  in  the  Bay,  but  the  interval  ended  on  station  17,  where  the  ship  maintained 
a  posititin  with  its  bow  into  the  wind,  which  was  from  the  north  at  16  to  24  knots.  While  the  ship  was 
on  this  station,  the  aerosol  scattering  coefficient  continued  to  drop  from  0  40  x  10  '*/m  to 
0.16  X  10^''/m,  and  the  relative  humidity  remained  constant  at  about  50  .  The  changes  in  the  si/e  and 
area  distributions  are  shown  in  Fig.  41,  where  curve  1  in  each  plot  is  the  average  of  two  distributions 
taken  outside  the  Bay.  Curve  2  is  the  si/.e  distribution  measured  at  the  entrance  to  the  Bay.  Curves  3 
and  4  are  averages  of  three  distributions  upon  the  ship's  arrival  at  and  departure  from  station  17  respec¬ 
tively. 
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If  the  450-m-level  air-mass  trajectory  for  this  interval  can  be  trusted  to  indicate  the  origin  of  the 
surface  air,  then  the  air  had  moved  rather  rapidly  (in  the  prior  12  hours)  from  upper  New  ii  ork  State 
down  the  east  coast  to  the  Chesapeake  Bay.  The  observed  change  in  the  si/e  distribution  for  larger  par¬ 
ticles  was  probably  more  related  to  the  changing  trajectory  of  the  air  mass  and  only  incident ly  related  tti 
the  entrance  of  the  ship  into  the  Chesapeake  Bay 

The  rapid  change  in  the  scattering  coefficient  makes  this  an  interesting  interval  tor  looking  at  the 
extinction  in  more  detail.  The  measured  scattering  coefficient  correspomling  to  surves  1  ihiough  4  of 
Fig.  41  were  1.85,  0.74,  0.40,  and  0.16  x  10  ''/m  respective!)  The  values  of  the  extinction  coefficient 
calculated  from  the  four  si/e  distributions  shown  in  Fig.  41  'or  0  55  yim  in  wavelength  .irc  I  60.  0  70, 
0.40,  and  0.16  x  10  '’/m  respectively  anil  <irc  in  excellent  agreement  with  the  measured  v,ilucs  id  the 
scattering  coefficient. 
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DISCUSSION 

Summary  of  the  Hygroscopicity  Results 

In  the  subsection  from  page  10  to  page  15  we  described  the  theory  and  experimental  method  used 
to  determine  the  value  of  the  particle  activation  parameter  (S®)  from  measurements  of  the  supersatura¬ 
tion  at  which  a  particle  of  known  dry  size  is  activated.  In  this  subsection  we  summarize  and  discuss  the 
results  of  these  measurements. 

Figure  42  shows  the  results  of  our  measurements  of  the  dry  size  (zq)  and  critical  supersaturation 
(SJ  of  atmospheric  aerosol  particles.  Each  data  point  represents  the  values  of  5;.  and  zq  at  the  center 
of  an  S^-ro  block,  such  as  is  shown  in  Fig.  7.  These  measurements  were  made  at  Washington,  D.C,, 
during  1978  to  1980  and  off  the  east  coast  of  the  United  States  during  the  cruise  of  the  USNS  ffayes  in 
April  1980.  Figure  42  also  shows  the  relationship  between  S,.  and  zq  for  sodium  chloride  particles 
(S°  =  1.3),  ammonium  sulfate  particles  (B°  ~  0.65),  and  particles  having  a  value  of  0.2.  The  S^-rg 
relationships  for  other  electrolytes  commonly  found  in  atmospheric  aerosol  particles,  such  as 
ammonium  nitiate,  sodium  nitrate,  and  sodium  sulfate,  lie  between  the  curves  for  sodium  chloride  and 
ammonium  sulfate.  Since  the  measured  values  of  critical  supersaturation  are  higher  than  those  calcu¬ 
lated  for  pure  particles  of  the  salts  commonly  present  in  atmospheric  particulates,  we  may  conclude  that 
a  significant  fraction  of  the  dry  mass  of  the  atmospheric  aerosol  particles  was  insoluble.  If  we  assume, 
for  instance,  that  the  soluble  material  in  atmospheric  aerosol  particles  is  ammonium  sulfate,  then  a 
value  of  of  0.2  implies  that  70%  of  the  mass  of  the  dry  particles  is  insoluble. 

Although  critical  supersaturation  was  measured  for  only  two  or  three  values  of  particle  radius  for 
any  given  aerosol  sample,  there  is  some  indication  in  the  results  presented  in  Fig.  42  that  the  slope  of 
the  Sf-Zo  relationship  is  not  as  steep  as  predicted  by  theory  (Eq.  (15))  for  the  case  when  particle  com¬ 
position  is  independent  of  particle  size.  One  possible  explanation  of  these  results  is  that  the  mass  frac¬ 
tion  of  soluble  material  in  the  particles  increased  with  decreasing  particle  size,  At  this  time,  however, 
such  a  conclusion  must  be  regarded  as  tentative.  Also,  a  number  of  approximations  were  made  in 
deriving  Eq.  (14).  In  particular,  it  was  assumed  that  at  the  critical  supersaturation  zVzq  is  large  com¬ 
pared  to  unity.  This  approximation  becomes  less  tenable  as  rg  and  *  assume  small  values.  However, 
the  error  introduced  by  the  approxim.ation  will  not  be  large  when  rg  >  0.02  jutn  and  e  >  0.1. 
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DRY  SIZE,  tq 

Mg  42  —  Sumnijrv  «i[  (he  hviiroscopiciu  mcasuremcnis  in  comparison 
\\[\h  calculated  '.alucs  for  pure  particles  (25'’'(  I 

Table  3  lists  the  values  of  fl®  computed  from  Eq.  (15)  and  the  measured  values  of  and  r^.  The 
value  of  fl®  was  found  to  range  from  0.11  to  0.55,  with  the  average  of  all  24  values  being  0.32.  Except 
for  the  data  of  21  April  the  indicated  particle-size  range  encompasses  two  or  three  size  ranges  transmit¬ 
ted  by  the  mobility  analyzer.  While  0.55  is  the  largest  value  of  we  have  measured  to  date,  Hanel 
and  Lehmann’s  (1981)  data  on  po  and  t)“  (Fig.  5)  for  aerosol  samples  yield  values  of  as  high  as  0  98 
for  particles  0.15  pm  to  5.0  pm  in  radius  which  were  90%  sea  salt.  The  number  of  measurements  of 
B°  which  we  have  made  is  not  sufficient  to  allow  us  to  draw  conclusions  regarding  differences  in  the 
value  of  B°  between  continental  and  maritime  air  for  particles  smaller  than  0.1  pm.  However,  if  the 
ocean  surface  is  not  a  significant  source  of  particles  in  this  size  range,  then  we  would  not  expect  large 
systematic  differences  in  fl®  between  continental  and  maritime  air  masses. 


Table  3  —  Values  of  B^  for 
atmospheric  aerosol  particles 


Location  of 
Measurements 

— 

Date 

Particle  Size 
Range  Measured 
(pm) 

Washington,  D.C. 

29  June  1978 

0.020-0.052 

0.23-0.40 

Washington,  D.C. 

26  July  1979 

0.043-0.079 

0.11-0.39 

Washington,  D  C. 

27  July  1979 

0.029-0.079 

0.19-0.27 

Washington,  D.C. 

30  July  1979 

0.029-0.079 

0.16-0.29 

Washington,  D.C. 

29  Aug.  1980 

0.020-0.035 

0.29-0.38 

Washington,  D  C. 

8  Sept.  1980 

0  020-0.035 

0.44-0.46 

200  km  off  the  Virginia  coast 

17-19  Apr.  1980 

0.015-0.047 

0.28-0.42 

Norfolk  harbor 

21  Apr.  1980 

0.032-0.047 

0.51 

Chesapeake  Bay 

22  Apr.  1980 

0.015-0.032 

0.40-0.55 

Figure  43  shows  the  dependence  of  particle  size  on  relative  humidity  for  =  0.1,  0.3,  and  0.55. 
These  curves  were  computed  from  Eqs.  (5),  (10),  (12),  and  (15)  using  B  =  0.8fl“  and  /)  =  20.2 
;uN/mm.  The  figure  shows  that  the  observed  range  in  gives  rise  to  large  differences  in  the  equili¬ 

brium  size  of  particles  at  high  relative  humidity.  The  figure  also  shows  the  effect  of  dry  particle  radius 
on  the  particle  growth  curves.  If  we  assume  simply  that  the  aerosol  extinction  coefficient  varies  as  the 
square  of  particle  radius,  then  a  factor-of-2  difference  in  particle  radius  means  a  factor-of-4  difference  in 
visibility. 


Relationship  Between  the  CCN  Supersaturation  Spectrum  and  the  Size  Distribution 

There  is  a  unique  relationship  between  the  dry  size  of  a  particle  and  its  critical  supersaturation  5^. 
This  relationship  can  be  formulated  in  terms  of  the  activation  parameter  and  is  given  by  Eq.  (15).  It 
is  plotted  for  various  values  (chemical  compositions)  in  Fig.  42.  By  use  of  the  expression  given  in 
Eq.  (15)  it  is  possible  to  convert  all  size  distributions  to  CCN  supersaturation  spectra.  The  supersatura¬ 
tion  spectrum  is  defined  as  the  cumulative  number  of  particles  activated  at  a  given  supersaturation; 
therefore  the  differential  size  distribution  must  first  be  integrated  to  obtain  the  cumulative  size  distribu¬ 
tion.  Curves  1,  2,  and  3  of  Fig.  44  show  supersaturation  spectra  calculated  for  a  size  distribution  taken 
on  19  April  at  0900  EDT  assuming  values  of  fl®  of  0.1,  0.2,  and  0.3,  respectively.  The  CCN  concentra¬ 
tions  taken  concurrently  are  plotted  as  individual  points  and  lie  below  those  calculated  from  the  size 
distribution.  The  degree  of  agreement  shown  between  calculated  CCN  concentrations  and  those  meas¬ 
ured  with  the  TGDCC  are  typical  of  those  obtained  during  the  cruise.  The  shape  of  the  two  curves  in 
the  overlap  region  usually  agreed  better  than  the  number  concentration.  Another  comparison  between 
measured  and  calculated  CCN  supersaturation  spectra  for  different  times  on  the  cruise  is  given  by 
Fitzgerald  and  Hoppel  (1982). 

The  measured  concentration  of  CCN  was  nearly  always  less  than  the  concentration  inferred  from 
the  size  distribution.  In  the  past  (Hoppel  and  Wojciechowski,  1976),  we  have  observed  that  the  CCN 
concentrations  obtained  from  the  video  recording  were  about  10  to  15"/a  below  that  obtained  from 
simultaneous  photographs.  The  photographic  method  has  a  smaller  minimum  detectable  size.  The 
TGDCC  data  presented  in  this  report  were  taken  with  the  videotape  system  Other  limitations  of  the 
TGDCC  also  tend  to  make  the  measured  CCN  concentrations  on  the  low  side  (Hoppel  and 
Wojciechowski,  1976). 
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Figure  45  shows  the  correlation  of  the  concentrations  of  CCN  active  at  0.7%  supersaiuraiion  with 
the  number  of  particles  in  channels  larger  than  0.0318  Mtn  obtained  from  the  size  distribution  taken 
closest  in  time  to  the  CCN  data.  Figure  45  shows  the  high  concentrations  of  CCN  observed  during  the 
cruise.  The  majority  of  the  concentrations  are  greater  than  2000/cm^  and  account  for  a  large  fraction  of 
the  total  aerosol  concentration  on  many  occasions.  The  reason  CCN  concentrations  are  often  a  large 
fraction  of  the  total  concentration  is  obvious  from  the  size  distributions  and  is  a  result  of  the  small 
number  of  particles  at  radii  below  0.04  /am  in  air  masses  that  have  been  over  water  for  a  number  of 
hours.  These  conditions  over  water  are  in  strong  contrast  to  conditions  usually  found  over  land.  Large 
CCN  concentrations  were  also  observed  on  a  1977  trans-Atlantic  cruise  as  the  ship  left  Norfolk  and 
moved  northeast  along  the  east  coast.  On  the  1977  cruise  the  CCN  count  remained  high  until  the  ship 
was  beyond  Nova  Scotia,  where  it  encountered  a  polar  air  mass  (Hoppel,  1979).  The  CCN  count  at 
that  point  dropped  to  levels  well  below  400/cm^. 
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It  is  generally  accepted  that  the  effects  of  water-vapor  depletion  in  the  TGDCC  due  to  growing 
droplets  is  small  for  concentrations  of  CCN  below  1000/cm'  As  the  concentration  increases,  vapor 
depletion  becomes  an  increasingly  serious  problem  (Chodes  et  al.,  1974,  and  others).  It  is  evident 
from  Fig.  45  that  the  measured  concentration  of  CCN  tends  '.o  be  smaller  than  the  number  of  particles 
larger  than  0.0318  nm  in  radius  when  the  CCN  concentration  is  large.  Since  the  direct  measurements 
of  would  indicate  that  particles  as  small  as  0.025  /am  should  be  activated  at  0.7%  supersaturation,  the 
CCN  concentration  should  be  at  least  as  large  as  the  concentration  of  particles  larger  than  0.0318  ^im. 
The  low  CCN  concentration  at  larger  concentration  is  probably  the  result  of  vapor  depletion  in  the 
TGDCC  causing  a  depressed  supersaturation.  If  we  attempt  to  determine  by  matching  the  supersa¬ 
turation  spectrum  inferred  from  the  size  distribution  with  that  from  the  TGDCC,  the  value  of  B“  would 
be  considerably  smaller  than  that  determined  from  the  more  direct  method  used  in  the  preceding  sub¬ 
section.  We  believe  the  more  direct  method  is  more  accurate,  because  it  does  not  depend  on  the  abso¬ 
lute  accuracy  of  the  measured  CCN  and  particle  concentrations  but  only  on  the  accuracy  of  the  size 
interval  transmitted  by  the  mobility  analyzer  and  the  accuracy  of  the  measurement  of  the  critical  super¬ 
saturation  in  the  TGDCC. 

Mechanisms  Which  Determine  the  Aerosol  Size  Distribution  in 
Air  Masses  Advecting  Off  the  East  Coast 

In  principle  it  should  be  possible,  starting  with  a  given  size  distribution  at  some  point  on  an  air- 
mass  trajectory,  to  predict  the  future  size  distribution  along  that  trajectory.  In  practice  this  is  generally 
not  possible,  because  one  cannot  accurately  quantify  the  sources,  sinks,  and  transformations  of  aerosol 
material,  and  the  limitations  of  models  to  accurately  predict  the  trajectory  and  meteorological  conditions 
along  that  trajectory.  The  very  simple  case  of  an  air  mass  advecting  directly  off  the  coast  in  the  absence 
of  precipitation  and  fronts  is  a  possible  exception.  In  this  '  tter  case  the  rapid  initial  decay  of  small  par¬ 
ticles  during  the  first  24  hours  can  be  explained  by  coagulation  of  particles  less  than  0.03  /am  with 
larger  particles  and  by  growth  of  small  particles  due  to  condensation  of  newly  formed  gas-phase  reaction 
products  of  low  volatility  on  their  surfaces.  An  example  of  the  decay  of  small  particles  as  an  air  mass 
moves  off  the  coast  is  shown  in  Fig.  31. 

Once  the  meteorological  processes  become  more  complex,  the  quantitative  interpretation  of  the 
observed  changes  in  the  size  distribution  becomes  difficult.  Some  of  the  important  but  difficult 
processes  are:  the  effect  of  a  capping  inversion  on  the  boundary  layer,  mixing  within  the  planetary 
boundary  layer,  the  cycling  of  air  through  clouds  which  subsequently  evaporate  (this  mechanism 
perhaps  having  an  important  effect  on  the  size  distribution  and  to  be  discussed  later),  removal  of  aero¬ 
sols  by  precipitation  scavenging,  formation  of  new  particles  by  gas-to-particle  conversion,  and  the  gen¬ 
eration  of  new  particles  at  the  sea  surface  during  high  seas  and  strong  wind.  Presently  there  are  few 
marine  aerosol  data  and  supporting  meteorological  data  available  to  sort  out  the  relative  contributions  of 
the  various  mechanisms. 

Our  observations  in  air  masses  advecting  off  the  east  coast  (in  the  absence  of  a  precipitation  cycle) 
are  consistent  for  the  most  part  with  what  would  be  expected  from  aerosol  theory: 

•  There  is  a  rapid  decay  of  small  particles  (<0.05  /am)  during  the  first  hours  as  the  air  mass 
moves  out  over  the  ocean  away  from  local  sources  of  particulates.  The  time  associated  with 
this  decay  is  such  that  the  particle  lifetime  increases  as  the  size  increases.  For  particles  0.05 
/am  in  radius  and  larger,  coagulative  decay  was  not  observed  over  the  interval  of  about  12 
hours. 

•  The  ab.sence  of  local  sources  makes  the  si/c  distribution  remarkably  stable,  often  over  the 
interval  of  hours  (after  the  initial  decay  of  small  particles)  Although  this  observation  is 
entirely  reasonable,  the  observed  stability  is  quite  remarkable  to  anyone  experienced  in 
measuring  urban  aerosols.  The  temporal  stability  of  the  size  distribution  measured  from  a 
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moving  platform  in  an  atmosphere  in  which  there  is  obviously  turbulent  mixing  implies  a 
degree  of  spatial  uniformity  as  well.  This  spatial  homogeneity  must  have  extended  at  least  10 
to  hundreds  of  kilometers  horizontally  and  10  to  hundreds  of  meters  vertically  on  some  occa¬ 
sions. 

•  Particles  in  the  range  0.08  to  0.5  nm  in  radius  were  less  variable  than  the  smaller  particles. 
This  region  of  the  size  distribution  has  aptly  been  referred  to  as  the  accumulation  mode, 
because  the  mass  lost  by  coagulation  from  smaller  sizes  is  transferred  to  this  region  of  the 
distribution  within  a  time  comparable  to  the  lifetime  of  particles  in  the  troposphere.  Even 
though  the  transfer  of  particles  to  the  accumulation  mode  greatly  reduces  the  concentration 
of  small  particles,  the  mass  transfer  is  small  and  produced  no  consistent  observable  effect  on 
the  size  range  0.08  to  0.5  /zm,  over  the  hours  when  the  decay  in  the  small-size  range  was 
greatest.  The  temporal  stability  of  the  accumulation  mode  also  implies  that  ocean  surface 
deposition  of  particles  in  this  size  range  is  small,  as  would  be  expected, 

The  decay  time  for  aerosols  smaller  than  0.05  /zm  strongly  depends  on  their  size,  as  is  evident 
from  Fig.  31,  and  it  makes  little  sense  to  discuss  a  residence  or  lifetime  of  atmospheric  aerosols  without 
specifying  the  size  range  of  interest.  For  example,  the  total  particulate  count  dropped  from  about 
9000/cm^  to  5500/cm^  between  curves  2  and  4  of  Fig.  31,  during  which  time  the  air  had  spent  7  hours 
over  water.  This  would  imply  a  rather  short  lifetime  for  aerosols.  This  loss  occurred  almost  entirely  in 
the  size  range  smaller  than  0.03  /zm.  There  was  no  observable  loss  of  mass  and  only  a  small  change  in 
the  measured  CCN  concentration.  For  particles  between  0.08  and  0.5  /zm  in  radius  the  lifetime 
depends  less  on  the  properties  of  the  aerosol  and  more  on  the  meteorological  conditions  to  which  the 
air  is  subjected;  for  example,  these  particles  are  removed  by  wet-scavenging  processes. 

Possible  Mechanisms  Which  Would  Produce  the  Double'Peaked  Size  Distributions 

On  several  occasions  size  distributions  occurred  with  prominent  double  peaks  in  the  small-particle 
size  range  (below  0.1  /zm)  and  were  quite  constant  for  a  number  of  hours.  Averages  during  two  of 
these  occasions  were  shown  as  curve  1  in  the  left-hand  plot  of  Fig.  38  and  curve  1  in  the  left-hand  plot 
of  Fig.  39;  each  individual  distribution  of  the  average  displayed  the  double  peak.  Figure  46  shows  a 
composite  of  five  of  the  eight  size  distributions  from  which  the  average  shown  as  curve  1  in  the  left- 
hand  plot  of  figure  38  was  calculated.  We  can  suggest  at  least  four  physical  mechanisms  which  might 
produce  a  double-peaked  size  distribution:  mixing  of  two  air  masses  with  very  different  aerosol  popula¬ 
tions,  heterogeneous  nucleation  and  growth  of  larger  particles  by  condensation  of  gas-phase  reaction 
products,  homogeneous  heteromolecular  nucleation  of  new  small  particles  in  air  initially  deficient  of 
particles  (the  peak  at  the  smaller  size  resulting  from  newly  formed  small  particles  and  the  peak  at  the 
larger  sizes  representing  the  remnant  of  the  size  distribution  in  an  aged  air  mass),  and  increase  in  mass 
of  CCN  which  have  been  through  a  cloud  condensation-evaporation  cycle. 

The  first  of  the  mechanisms  which  might  produce  a  double  peak  is  the  mixing  of  two  air  masses 
with  very  different  size  distributions,  ft  is  certainly  possible  to  invent  two  size  distributions  which, 
when  added  together,  will  result  in  a  size  distribution  such  as  those  shown  in  Fig.  46.  However,  the 
two  size  distributions  (assuming  each  is  singly  peaked)  must  be  much  different  from  what  we  have  ever 
observed.  To  reproduce  the  strong  minimum,  the  size  distribution  exhibiting  the  peak  at  the  smaller 
size  must  decrease  with  increasing  size  much  more  rapidly  than  has  ever  been  observed.  Therefore,  it 
is  highly  unlikely  that  the  size  distribution  of  Fig.  46  was  formed  by  the  mixing  of  air  masses. 

The  second  mechanism  which  might  produce  a  double  peak  is  the  continual  formation  of  new 
gaseous  species  in  the  atmosphere  by  gas-phase  reactions  (including  photochemical  reactions).  Some  of 
these  reaction  products  are  of  low  volatility  and  will  diffuse  to  and  condense  on  existing  particles.  If 
the  particles  are  larger  than  a  critical  size,  they  will  be  nucleated  and  grow  (details  of  the  mechanism 
being  given  in  Hoppel,  1975).  This  causes  a  separation  in  the  size  distribution  at  the  critical  radius  and 
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Fig  46  —  C  omposilc  of  five  consccuti\c  si/e  distribu¬ 
tions  over  ii  3-hour  interval  which  all  exhibited  a  dou¬ 
ble  peak 


RADIUS  turn! 


could  lead  to  the  double-peaked  size  distribution.  Use  of  estimated  gas-phase  reaction  rates  for  the 
marine  environment  to  quantitatively  estimate  the  growth  rate  of  particles  with  radii  from  0.03  to  0.06 
34m  by  this  mechanism,  however,  yields  growth  times  which  are  much  longer  than  the  lifetimes  of 
atmospheric  aerosols.  It  is  therefore  doubtful  that  this  mechanism  produced  the  size  distributions 
shown  in  Fig.  46. 


The  third  mechanism  which  might  produce  a  double  peak  is  the  production  of  new  particles  at 
small  sizes  (<0.005  /xm).  They  are  produced  when  the  surface  area  of  the  existing  particles  is 
insufficient  to  absorb  the  newly  formed  involatile  products  and  the  concentration  of  these  involatile  pro¬ 
ducts  increases  to  the  level  required  for  spontaneous  nucleation  of  new  particles.  These  newly  formed 
particles  grow  by  condensation  of  additional  gas-phase  reaction  products  and  by  coagulation  with  other 
particles.  This  growth  becomes  slow  as  the  particle  radius  increases  above  0.01  4m.  The  newly  formed 
particles  produce  a  new  peak  at  a  smaller  radius,  in  addition  to  the  residual  peak,  and  could  certainly 
produce  a  double-peaked  size  distribution.  This  mechanism  has  been  suggested  by  others  and  possibly 
has  been  observed  by  Haaf  and  Jaenicke  (1980)  on  a  mountain  site  (Schauinsland  in  the  Black  Forest), 
where  the  increase  in  small  particles  in  the  morning  was  strongly  correlated  with  solar  radiation.  Haaf 
and  Jaenicke  found  a  weak  suggestion  in  their  data  of  small  particle  formation  in  the  eastern  North 
Atlantic  on  sunny  days.  We  looked  in  vain  for  a  correlation  between  small  particles  and  solar  radiation 
throughout  the  /^ayes  cruise  and  could  find  none.  If  there  were  solar  influenees  on  the  concentration  of 
small  particles,  they  were  not  immediate  and  were  masked  by  the  more  prominent  meteorological 
effects.  Furthermore  the  total  aerosol  load  off  the  east  coast  of  the  United  States  during  our  cruise  was 
considerably  greater  than  that  usually  encountered  in  the  middle  and  eastern  Atlantic.  Thus  the  surface 
area  available  to  absorb  newly  formed  gas-phase-reaction  products  is  larger  than  in  more  remote  regions 
of  the  Atlantic.  In  regions  where  the  aerosol  load  is  less,  the  probability  of  new  particle  formation  is 
greater.  The  double-peaked  size  distribution  shown  in  Fig.  46  exhibited  few  particles  with  a  radius  less 
than  0.01  /xm.  This  would  indicate  that  if  nucleation  of  new  particles  occurred,  then  the  mechanism 
acted  earlier  and  new  particle  formation  had  either  ceased  or  was  weak  at  the  time  the  measurements 
were  made.  These  measurements  were  made  on  a  clear  sunny  morning. 

The  fourth  mechanism  which  might  produce  a  double  peak  is  the  enhancement  of  involatilc 
material  in  a  droplet  during  formation  and  reevaporation  of  cloud  droplets.  Pruppachcr  and  KIctt 
(1978)  estimate  that  globally  the  residence  time  of  water  in  the  troposphere  is  about  9  days  before 
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removal  by  precipitation  and  that  there  are  about  ten  condensation-evaporation  cycles  before  a  precipi¬ 
tating  system  is  tormed.  This  would  indicate  that  moisture  goes  through  a  cloud  formation-and- 
evaporation  cycle  about  once  a  day.  In  a  well-mixed  lower  troposphere  this  would  imply  that  air  (with 
aerosols)  is  also  filtered  through  clouds  at  about  the  same  freciuency.  When  a  cloud  forms,  aerosols 
larger  than  a  critical  size  will  be  activated  and  grow  to  cloud  drop  sizes.  Ty  picall)  the  change  in  size  will 
be  from  0.05  /am  to  5.0  /am,  which  is  an  increase  in  mass  of  the  order  of  10*’.  The  increased  size  of 
cloud  droplets  makes  them  effective  surfaces  to  which  smaller  aerosol  particles  can  diffuse  and  effective 
centers  for  absorbtion  of  various  trace  gases.  Some  of  these  dissolved  gases  react  with  other  dissolved 
gases  and  water  to  form  new  substances  which  will  remain  as  additional  residue  when  the  cloud  droplet 
reevaporates.  The  process  of  this  type  most  widely  studied  is  the  conversion  of  SO2  to  sulfate  in  the 
cloud  droplet,  by  Easter  and  Hobbs  (1974)  and  Hegg  and  Hobbs  (1979).  If  only  one  foreign  molecule 
in  lO*"  condensed  water  molecules  is  retained  as  new  nuclear  material,  then  upon  reevaporation  of  the 
droplet,  the  mass  of  the  residue  will  be  increased  by  as  much  as  a  factor  of  2  and  the  radius  will  be 
increased  by  the  cube  root  of  2  (being  less  for  larger  CCN). 

To  illustrate  how  this  process  might  work  to  form  a  double-peaked  size  distribution,  a  typical  size 
distribution  shown  as  curve  1  in  Fig.  47  has  been  subjected  to  a  simulated  1-hour  cloud  cycle.  The 
maximum  supersaturation  in  the  cloud  is  assumed  to  be  about  0.3%,  so  that  all  particles  greater  than 
about  0.04  /zm  are  activated  to  form  a  cloud  with  a  droplet  concentration  of  500/cm’  and  a  liquid  water 
content  of  0.3  g/m’.  Nuclei  smaller  than  0.04  /zm  are  not  activated  and  are  allowed  to  diffuse  to  the 
cloud  droplets  according  to  the  following  formula  by  Twomey  (1977): 

(2J) 

nQ(r) 

where  n(r)  is  the  number  of  particles  of  radius  r(r  <  0.04 /am)  that  diffuse  to  the  cloud  droplets, 
no(r)  is  the  initial  number  of  particles  of  radius  r,  IE  is  the  liquid  water  content  of  the  cloud,  /V  is  the 
cloud  droplet  concentration,  D(r)  is  the  diffusion  coefficient  for  particles  of  radius  r,  and  i  is  the  time. 
Upon  reevaporalion,  particles  initially  bigger  than  0.04  /am  in  radius  were  assumed  to  have  arbitrarily 
increased  their  mass  by  a  factor  of  2,  due  to  the  in-cloud  processes  mentioned.  Curve  2  of  Fig.  47  is 
the  size  distribution  obtained  from  curve  1  when  unactivated  particles  are  lost  to  cloud  droplets  by 
diffusion  according  to  Eq.  (21)  and  the  mass  of  cloud  nuclei  has  been  increased  by  a  factor  of  2  during 
the  condensation/evaporation  cycle.  The  depth  of  the  minimum  between  the  double  peak  of  curve  2 
was  drawn  in  arbitrarily.  Ideally  the  minimum  would  be  even  less,  but  any  mixing  with  air  which  had 
not  been  through  a  cloud  cycle  has  a  strong  tendency  to  destroy  the  minimum.  If  this  mechanism  is 
responsible  for  the  double-peaked  size  distribution  of  Fig.  46,  then  the  strong  minimum  observed 
implies  that  the  air  had  been  through  a  cloud  cycle  at  least  once  and  was  not  mixed  with  air  which  had 
not  been  through  a  cloud  cycle. 

The  size  distributions  of  Fig.  46  were  obtained  under  clear  conditions  just  behind  a  cool  front 
However,  the  NOAA  back-trajectories  indicated  that  the  air  had  originated  in  a  cloudy  overcast  region 
just  behind  the  front  within  the  prior  12  hours.  The  reliability  of  the  trajectories  under  these  conditions 
is  unknown.  It  is  interesting  that  Haaf  and  Jaenicke  (1980)  observed  on  Schauinsland  (at  1230  m)  on 
"days  with  complete  cloud  cover"  size  distributions  somewhat  similar  to  some  of  our  double-peaked  dis¬ 
tributions. 


Of  the  four  mechanisms  discussed  which  may  produce  double-peaked  distributions,  the  lalter  two 
are  probably  both  important  at  times  in  the  marine  atmosphere  However,  our  discussion  is  specula¬ 
tive,  and  more  aero.sol  measurements  and  supporting  meteorological  measurements  are  needed  to 
resolve  specific  mechanisms. 
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Summary  of  Electromagnetic  Extinction  Results 


The  exlinciion  coefficienis  al  0.55  in  wavelength  for  those  size  distributions  from  the  1980 
USNS  Hayes  cruise  that  are  given  in  the  figures  of  this  report  (excluding  composites)  are  tabulated  in 
Table  4.  The  extinction  is  computed  by  use  of  approximation  formulas  for  the  Mie  theory  extinction 
efficiency  factor  Q,.  (fitzgerald,  1979),  analytic  approximations  to  the  particle  distributions  in  the  form 
of  the  sum  of  three  log-normal  distributions,  index-of-refraction  data  frr)m  Paltridge  and  Platt  (1976) 
and  Shetlle  and  Fenn  (1979),  and  the  relative  amounts  of  water  and  dry  material  in  the  particles  as  a 
function  of  relative  humidity,  as  given  by  Fitzgerald  (1978).  The  details  of  the  compulations  were  dis¬ 
cussed  in  the  subsection  beginning  on  page  15. 


The  agreement  between  measured  and  computed  extinction  as  given  in  Table  4  is  good  for  the 
intervals  when  the  ship  was  under  high  pressure  off  the  cast  coast  (Fig.  .18.  curve  1).  and  for  the  inter¬ 
val  of  rapid  change  in  extinction  coefficients  just  before  the  return  to  Chesapeake  Bay  (Fig.  41).  The 
agreement  was  also  good  for  the  interval  when  the  distributions  were  stable  and  double-peaked  (Fig. 
39,  curve  1).  Further,  although  the  computed  values  of  extinction  for  before  and  after  passage  of  a 
cold  front  on  21  April  (Fig.  40)  were  lower  than  the  measured  extinctions,  the  ratio  of  the  extinction 
after  passage  of  the  front  to  that  before  passage  is  nearly  the  same  for  measured  and  computed  extinc¬ 
tion. 


The  percentage  of  the  computed  extinction  at  0.55  /um  in  wavelength  due  to  particles  in  selected 
size  ranges  shown  in  Table  4  is  plotted  in  Fig.  48.  There  arc  wide  variations  in  the  contributions  of  the 
individual  size  ranges  to  the  total  extinction,  although  the  contribution  from  particles  with  a  radius  less 
than  0.1  /zm  is  small.  However,  the  variation  during  unchanging  weather  conditions  is  much  less. 
Only  in  the  rapidly  changing  air  masses  after  the  ship  left  Chesapeake  Bay  (Fig.  31)  and  in  the  changing 
air  masses  associated  with  passage  of  a  warm  front  (Fig.  34)  were  there  appreciable  changes  in  the  rela¬ 
tive  contribution  of  different  size  ranges,  as  can  be  seen  from  Table  4. 
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Table  4  —  Extinction  coefficients  at  0  55  in  wavelength 
for  particle  distributions  shown  in  this  report 


Figure-Curve 

Extinction 
(10  Vm) 

Rll' 

(%) 

Extitiction  (%)  in  Particle-Radius  Ranges' 

0.04  to 
0.1  /u  m 

0.1  to 
0.3  /am 

0.3  to 
0.6  /am 

Measured* 

Computed^ 

0.9  /am 

3  /XI 

31-1 

0.36 

0.28 

42  ^ 

7 

75 

7 

1  5 

■  6 

31-2 

0.55 

0.51 

65 

6 

38 

23 

16 

19 

31-3 

0.46 

0.33 

65 

7 

36 

17 

18 

31-4 

0.55 

0.39 

52 

8 

55 

9 

9 

18 

33-1 

1.76 

1.10 

86 

2 

30 

26 

23 

18 

33-2 

1.51 

1.54 

80 

2 

26 

33 

23 

16 

33-3 

1.01 

1.23 

80 

1 

18 

28 

30 

23 

33-4(34-1) 

0.80 

0.73 

72 

1 

23 

30 

24 

34-2 

— 

95 

1 

14 

31 

34 

34-3 

- 

1.69 

95 

4 

21 

33 

42 

38-1 

0.35 

64 

4 

47 

25 

12 

II 

38-1 

0.32 

6 

49 

31 

8 

6 

38-1 

0.19 

45 

5 

61 

23 

5 

5 

39-1 

0.31 

38 

7 

68 

18 

3 

3 

39-2 

0.64 

3 

79 

14 

2 

2 

0.46 

84 

1 

67 

22 

5 

4 

0.93 

75 

3 

23 

8 

6 

41-1 

1.8 

74 

2 

39 

35 

15 

9 

41-2 

4 

47 

34 

8 

6 

41-3 

0.36 

55 

8 

53 

29 

4 

6 

41-4 

0.16 

58 

5 

54 

29 

5 

7 

'The  measured  extinction  is  really  the  scattering  coefficient  measured  with  an  integrating  nephelomcter 

^he  total  extinction  and  extinction  due  to  particles  in  selected  size  ranges  were  computed  by  the  method  outlined  in  the  subs 
tion  beginning  on  page  IS 

‘The  relative  humidity  was  determined  from  dew  point  and  temperature 


l  ig  48  —  Contribution  ol  sarious  si/c  ranges  to  the  total 
aerosol  extinction  .it  II  S'  pm  in  w.isclengih 
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Particles  0.1  to  0.3  /um  in  radius  contributed  the  most  to  the  esiinciion  in  16  of  the  21  distribu¬ 
tions.  The  least  variability  in  percent  contribution  to  total  extinction  was  observed  in  the  range  0.3  to 
0.6  juni.  This  is  consistent  with  I'igs.  27.  29,  and  30.  which  show  that  the  variations  in  the  size  distri¬ 
bution  at  intermediate  sizes  are  much  less  than  the  v  iriations  at  smaller  and  larger  sizes.  In  only  one 
distribution  did  particles  0.3  to  0.6  /um  in  radius  make  the  largest  contribution  to  total  extinction. 

The  contribution  to  total  extinction  at  O.-S.^  ^m  from  particles  greater  than  0.6  /am  in  radius  was 
also  highly  variable  from  case  to  case,  with  a  minimum  of  4%  and  a  maximum  cf  75%  of  the  extinction. 
The  contribution  from  particles  0.6  to  0.9  /am  in  radius  and  0.9  to  3  /am  in  radius  were  usually  nearly 
equal  and  tracked  each  other  in  their  contribution  to  the  total  extinction.  This  can  be  attributed  to  the 
origin  of  many  of  the  particles  in  these  size  ranges  being  from  wind  and  wave  action.  The  two  distribu¬ 
tions  in  which  particles  0.9  to  3  /am  in  radius  made  the  largest  contribution  to  total  extinction  were 
composites  of  distributions  in  high  winds,  as  was  discussed  with  regard  to  the  distributions  shown  in 
Fig.  34. 

Waggoner  et  al.  (1981)  have  reported  a  good  correlation  between  fine-particle  mass  and  visibility 
for  rural  and  urban  air.  To  test  this  correlation  for  the  //oves  cruise,  the  mass  was  computed  from  the 
log-normal  fit  to  the  data  points.  A  particle  density  of  2.0  g/cm'  and  a  large-radius  cutoff  of  1.25  /xm 
were  selected  for  comparison  with  the  results  of  Waggoner  et  al.  (1981). 

The  correlation  between  the  fine-particle  mass  and  extinction  at  0.55  /xm  is  shown  in  Fig.  49, 
where  the  extinction  as  measured  with  a  nephelometer  is  plotted  against  the  mass  of  particles  of  radius 
between  0.05  and  1.25  /xm  calculated  from  the  size  distributions.  The  correlation  coefficient  is  0.89, 
and  the  ratio  of  mass  concentration  to  extinction  is  0.60  g/m’,  which  is  within  the  range  of  0.2  to  0,7 
g/m^  indicated  by  the  data  of  Charlson  et  al.  (1968). 


Fig.  49  —  Correlation  between  fine-particle 
mass  (r  <  1.25  /xm)  and  extinction  (scattering 
coefficient)  at  0.55  /xm 


0  50  100 
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The  correlation  coefficient  between  the  extinction  computed  from  the  size  distribution  and  the 
mass  of  the  particles  (0.05  <  r  <  1.25  /xm)  is  0.98,  whereas  the  correlation  between  the  calculated  and 
measured  extinction  at  0.55  /xm  in  wavelength  is  0.91.  The  mass  is  relatively  independent  of  the  lower 
radius  cutoff,  with  a  change  in  the  lower  radius  cutoff  of  from  0.01  to  0.1  /xm,  resulting  in  a  few  per¬ 
cent  change  in  mass.  The  mass  is  much  more  dependent  on  the  large-radius  cutoff,  with  the  increase  in 
mass  resulting  from  an  extension  of  the  large-radius  cutoff  to  3  /xm  ranging  from  less  than  20"/(i  to  well 
over  100%. 
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Waggoner  el  al.  (1981)  also  present  a  historic  perspective  to  the  optical  characteristics  of  atmos¬ 
pheric  aerosols  before  discussing  the  results  of  their  measurements  of  aerosols  in  rural  and  urban  air 
The  correlations  between  fine-particle  mass  and  e.xtinction  for  their  measurements  were  greater  than 
0,95.  and  their  mass-to-exiinction  ratio  was  about  0.3  g/m’  for  the  continental  aerosols  they  measured 


CONCLUSIONS 

In  this  report  we  have  described  the  instruments  and  methods  developed  at  NRL  to  measure  sub¬ 
micrometer  particles  and  have  presented  the  results  of  size  distributions  and  solubility  measurements 
along  the  east  coast  of  the  United  Stales  both  at  a  shore  site  (Wallops  Island)  and  from  a  ship  within 
300  km  of  the  coast.  Some  measurements  made  at  Washington,  D  C.,  are  included  to  contrast  submi¬ 
crometer  aerosol  properties  in  urban  air  with  those  in  aged  continental  air  off  the  east  coast. 

Earlier  measurements  of  the  marine  size  distribution  were  of  such  poor  quality  that  it  was  impos¬ 
sible  to  determine  if  different  distributions  were  real  or  an  artifact  of  the  measurement  technique  (Jaen- 
icke,  1978).  The  size  distributions  presented  in  this  report  are  sufficiently  accurate  and  numerous  to 
assure  that  the  observed  changes  were  real  and  that  the  data  can  be  used  with  confidence  in  the  search 
for  physical  mechanisms  causing  the  changes. 

Some  of  the  observed  changes  can  be  understood  in  general  qualitative  terms: 

•  The  rather  rapid  decay  of  small  particles  as  an  air  mass  leaves  land  and  advects  over  the 
ocean  is  due  mainly  to  coagulation  of  small  particles  with  particles  greater  than  0.06  ^m  in 
radius  and  indicates  little  to  no  production  of  new  submicrometer  particles; 

•  The  observed  termporal  stability  of  the  size  distribution  over  many  hours  would  be  expected 
in  air  masses  in  which  the  initial  decay  of  small  particles  had  already  occurred  and  no  wet¬ 
scavenging  of  particles  is  occurring; 

•  The  decrease  in  particles  in  the  range  0.08  to  1.0  ^tm  in  an  air  mass  which  appeared  to  go 
through  a  precipitation  cycle  can  qualitatively  be  explained  in  terms  of  wet-scavenging  of  par¬ 
ticles; 

•  An  increase  in  particles  larger  than  0.5  /zm  during  heavy  seas  and  high  winds  is  undoublcdlv 
caused  by  surface-generated  sea-salt  particles. 

Many  observed  changes  in  the  size  distribution  are  difficult  to  explain  even  qualitatively  Some  of 
the  difficulty  in  understanding  the  changes  may  originate  with  our  inability  to  determine  the  history  of 
the  air  encountered  and  the  role  which  condensation  and  reevaporation  of  cloud  droplets  play.  The 
existence  of  two  distinct  peaks  on  several  occasions  is  an  interesting  example  of  a  phenomenon  which 
we  were  unable  to  explain  adequately. 

Efforts  to  quantify  or  model  mechanisms  believed  to  be  responsible  for  the  changes  have  not  yet 
been  attempted  for  specific  cases  given  in  this  report.  In  principle  it  should  be  possible  to  start  with  a 
size  distribution  at  one  point  on  the  air  mass  trajectory  and  predict  the  size  distribution  at  a  later  time 
In  practice  this  task  is  extremely  difficult,  because  a  prior  size  distribution  on  the  same  trajectorv  is  not 
usually  known  and  the  trajectories  and  the  meteoroligical  history  along  a  trajectory  cannot  be  predicted 
accurately.  Nevertheless  a  comparison  of  model  predictions  including  various  known  mechanisms 
applied  to  specific  cases  given  by  our  data  could  help  sort  out  the  relative  importance  of  mechanisms. 
Such  studies  are  planned. 
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Siiuc  the  measurenicnts  given  in  this  report  were  all  taken  within  300  km  of  the  coast,  the  air 
masses  en.-ouniered  were-  not  suflieiently  aged  that  the  aerosol  k>ad  could  be  considered  to  represent 
"marine"  or  "background"  levels.  "Marine"  aerosol  and/or  "background"  aerosol  are  poorly  defined 
descriptions  olten  used  to  indicate  low  submicrometer  concentrations  found  at  locations  remote  from 
land  and  thus  from  anthropogenic  sources  of  pollution.  Measurements  in  air  masses  farther  removed 
from  land  are  necessary  to  observe  the  change  from  the  aged-continental  size  distributions  found  near 
the  east  coast  of  the  United  States  to  background  levels  believed  to  exist  in  more  remote  regions 
Measurements  farther  from  land  were  made  on  a  cruise  during  March  and  April  of  1983  and  will  be  an 
important  addition  to  the  measurements  reported  here. 

Measured  values  of  the  aerosol  solubility  (activation)  parameter  B"  ranged  from  0.1  to  0.55.  The 
effect  of  the  value  of  5''  on  the  growth  of  particles  as  a  function  of  relative  humidity  was  shown  in  Fig. 
43.  A  particle  having  5**  =  0.55  and  a  dry  radius  of  0.05  /am  will  swell  to  twice  its  dry  size  at  95  '/"  rela¬ 
tive  humidity,  whereas  a  particle  with  fi"  =  0.1  will  gtow  to  only  1.3  times  its  dry  size  at  95  "/  relative 
humidity.  If  it  is  assumed  that  the  soluble  component  of  atmospheric  particles  is  ammonium  sulfate 
and  that  the  particles  have  a  dry  density  of  1.8  g/cm\  then  a  B*'  value  of  0.1  implies  that  15  v  of  the 
dry  mass  of  the  particle  is  soluble  and  a  S"  value  of  0.55  implies  that  80"/'  of  the  mass  is  soluble. 

Electromagnetic  extinction  in  the  wavelength  range  0.3  to  12  /am  was  calculated  from  the  mea.s- 
ured  size  distributions  using  approximation  formulas  for  the  Mie  extinction  efficiency  factor.  The 
agreement  between  the  calculated  extinction  at  0.55  /am  and  the  measured  scattering  coefficient  is  good 
Particles  in  the  size  range  0.1  to  0.3  /am  made  the  largest  contribution  to  the  extinction  on  nearly  all 
occasions.  A  high  correlation  was  found  between  fine-particle  mass  (calculated  from  the  measured  size 
distribution  using  a  particle  density  of  2.0  g/cm'  and  a  large-radius  cutoff  of  1.25  /zm)  and  the  meas¬ 
ured  scattering  coefficient.  The  mean  value  of  the  ratio  of  mass  concentration  to  the  scattering 
coefficient  at  0.55  /am  was  0.60  g/m^ 
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